GLOBAL STRONG SOLUTION TO THE DENSITY-DEPENDENT 
INCOMPRESSIBLE VISCOELASTIC FLUIDS 
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Abstract. The existence and uniqueness of the global strong solution with small ini- 
tial data to the three-dimensional density-dependent incompressible viscoelastic fluids is 
established. The local existence and uniqueness of the global strong solution with small 
initial data to the three-dimensional compressible viscoelastic fluids is also obtained. A 
new method is developed to estimate the solution with weak regularity. Moreover, as a 
byproduct, we show the global existence and uniqueness of strong solution to the density- 
dependent incompressible Navier-Stokes equations using a different technique from [8]. 
All the results apply to the two-dimensional case. 
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1. Introduction 

Elastic solids and viscous fluids are two extremes of material behavior. Viscoelastic 
fluids show intermediate behavior with some remarkable phenomena due to their "elastic" 
nature. These fluids exhibit a combination of both fluid and solid characteristics, keep 
memory of their past deformations, and their behaviour is a function of these old deforma- 
tions. Viscoelastic fluids have a wide range of applications and hence have received a great 
deal of interest. Examples and applications of viscoelastic fluids include from oil, liquid 
polymers, mucus, liquid soap, toothpaste, clay, ceramics, gels, some types of suspensions, 
to bioactive fluids, coatings and drug delivery systems for controlled drug release, scaffolds 
for tissue engineering, and viscoelastic blood flow flow past valves; see [T0l[13l[35] for more 
applications. For the viscoelastic materials, the competition between the kinetic energy 
and the internal elastic energy through the special transport properties of their respective 
internal elastic variables makes the materials more untractable in understanding their be- 
havior, since any distortion of microstructures, patterns or configurations in the dynamical 
flow will involve the deformation tensor. For classical simple fluids, the internal energy can 
be determined solely by the determinant of the deformation tensor; however, the internal 
energy of complex fluids carries all the information of the deformation tensor. The inter- 
action between the microscopic elastic properties and the macroscopic fluid motions leads 
to the rich and complicated rheological phenomena in viscoelastic fluids, and also causes 
formidable analytic and numerical challenges in mathematical analysis. The equations of 
the density-dependent incompressible viscoelastic fluids of Oldroyd type ([27||28j) in three 



spatial dimensions take the following form |12|. 1231 [29] : 

pt + div(pu) = 0, (1.1a) 

(pu) t + div (pu (g> u) - pAu + VP(p) = div(pFF T ), (1.1b) 

F t + u-VF = VuF, (1.1c) 

divu = 0, (l.ld) 



where p stands for the density, u € M 3 the velocity, and F 6 M 3x3 (the set of 3 x 3 matrices) 
the deformation gradient. The viscosity coefficient p > is a constant. The increasing 
convex function P(p) = Ap 1 is the pressure, where 7 > 1 and A > are constant. Without 
loss of generality, we set A = 1 in this paper. The symbol (g> denotes the Kronecker tensor 
product and F T means the transpose matrix of F. As usual we call equation (jl.lap the 
continuity equation. For system (jl.ip . the corresponding elastic energy is chosen to be the 
special form of the Hookean linear elasticity: 

^(F) = ^|F| 2 , 

which, however, does not reduce the essential difficulties for analysis. The methods and 
results of this paper can be applied to more general cases. 

In this paper, we consider equations (jl.ip subject to the initial condition: 

(p,u,E)\ t=0 = (p Q (x),u (x),E (x)), xeR 3 , (1.2) 

and we are interested in the global existence and uniqueness of strong solution to the 
initial- value problem (|l.l[) - (|1.2p near its equilibrium state in the three dimensional space 
IR 3 . Here the equilibrium state of the system (jl.ip is defined as: p is a positive constant 
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(for simplicity, p = 1), u = 0, and F = I (the identity matrix in M 3x3 ). We introduce a 
new unknown variable E by setting 



By a strong solution, we mean a triple (p, u, F) satisfying (II. ip almost everywhere with 
initial condition (|1.2|) . in particular, (p, u, F)(-, t) 6 VF 2,9 , g 6 (3,6] for almost all £ > in 
this paper. 

When density p is a constant, system (jl.ip governs the homogeneous incompressible 
viscoelastic fluids, and there exist rich results in the literature for the global existence of 
classical solutions (namely in if 3 or other functional spaces with much higher regularity); 
see [5j [EJ [TH [HI [HI [20], [23] and the references therein. When density p is not a 
constant, the question related to existence becomes much more complicated and not much 
has been done. In [17] the authors considered the global existence of classical solutions 
in H s of small perturbation near its equilibrium for the compressible viscoelastic fluids 
without the pressure term. One of the main difficulties in proving the global existence 
is the lacking of the dissipative estimate for the deformation gradient. To overcome this 
difficulty, the authors in [16] introduced an auxiliary function to obtain the dissipative 
estimate, while the authors in [18] directly deal with the quantities such as Au + divF. 
Those methods can provide them with some good estimates, partly because of their high 
regularity of (u, F). However, in this paper, we deal with the strong solution with much 
less regularity in W 2,q , q 6 (3, 6], hence those methods do not apply. Thus, we need a new 
method to overcome this obstacle, and we find that a combination between the velocity 
and the convolution of the divergence of the deformation gradient with the fundamental 
solution of Laplace operator will develop some good dissipative estimates required for the 
global existence. The global existence of strong solution to the initial-value problem (jl.ip - 
(jl.2p is established based on the local existence and global uniform estimates. The local 
existence is obtained using a fixed point theorem without incompressible condition (jl.ldp . 
that is, the local existence holds for both the incompressible viscoelastic fluids (|l.la|) - 
(jl.ldp and the compressible viscoelastic fluids (jl.laj) - (|l.lcj) . The global existence and 
uniqueness of strong solution also holds for the density-dependent incompressible Navier- 
Stokes equations when the deformation gradient does not appear, which, as a byproduct, 
gives a similar result to [8] but through a different technique. 



The viscoelastic fluid system (II. ID can be regarded as a combination of the inhomo- 
geneous incompressible Navier-Stokes equation with the source term div(pFF T ) and the 
equation (ll.lcl) . For the global existence of classical solutions with small perturbation 
near an equilibrium for the compressible Navier-Stokes equations, we refer the reader to 
[24\ [25| [26| [30] and the references cited therein. We remark that, for the nonlinear invis- 
cid elastic systems, the existence of solutions was established by Sideris-Thomases in [33] 
under the null condition; see also [31] for a related discussion. 

The existence of global weak solutions with large initial data of (II. ip is still an out- 
standing open question. In this direction for the homogeneous incompressible viscoelastic 
fluids, when the contribution of the strain rate (symmetric part of Vu) in the constitutive 
equation is neglected, Lions-Masmoudi in [22] proved the global existence of weak solu- 
tions with large initial data for the Oldroyd model. Also Lin-Liu-Zhang in [19] proved the 



F = I + E. 



Then, (fLlc|) reads 



E t + u • VE = VuE + Vu. 



(1.3) 
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existence of global weak solutions with large initial data for the incompressible viscoelastic 
fluids when the velocity satisfies the Lipschitz condition. When dealing with the global 
existence of weak solutions of the viscoelastic fluid system (11. lh with large data, the rapid 
oscillation of the density and the non-compatibility between the quadratic form and the 
weak convergence are two of the major difficulties. 

The rest of the paper is organized as follows. In Section 2, we recall briefly the density- 
dependent incompressible viscoelastic fluids from some basic mechanics and conservation 
laws. In Section 3, we state our main results, including the local and global existence and 
uniqueness of the strong solution to the equations of the viscoelastic fluids, as well as to the 
incompressible Navier-Stokes equations with small data. In Section 4, we prove the local 
existence via a fixed-point theorem. In Section 5, we prove the uniqueness of the solution 
obtained in Section 4. In Section 6, we establish some global a priori estimates, especially 
on the dissipation of the deformation gradient and gradient of the density. In Section 7, 
we first prove some energy estimates uniform in time and some refined estimates on the 
density and the deformation gradient, and then give the proof of the global existence. 



2. Background of Mechanics for Viscoelastic Fluids 

To provide a better understanding of system (jl.ip . we recall briefly some background 
of viscoelastic fluids from mechanics in this section. 

First, we discuss the deformation gradient F. The dynamics of a velocity field u(x,t) 
in mechanics can be described by the flow map or particle trajectory x(t, X), which is a 
time dependent family of orientation preserving diffeomorphisms defined by: 

f f t x(t,X)=u(t,x(t,X)), 
\ x(0,X)=X, 

where the material point X (Lagrangian coordinate) is deformed to the spatial position 
x(t, X) , the reference (Eulerian) coordinate at time t. The deformation gradient F is 
defined as 

F(t,X) = ^(t,X), 

which describes the change of configuration, amplification or pattern during the dynamical 
process, and satisfies the following equation by changing the order of differentiation: 

9F(t,X) _ du(t,x(t,X)) 

dt ~ dx [ ] 

In the Eulerian coordinate, the corresponding deformation gradient F(t, x) is defined as 

F(t,x(t,X))=F(t,X). 
Equation (|2.2[) . combined with the chain rule and (|2.1|) . gives 

d t F(t, x(t, X))+u- VF(t, x(t, X)) = d t F(t, x(t, X)) + g F(^(^)) . ^Ml 

_ dF(t,X) _ du(t,x(t,X)) _ du(t,x(t,X)) dx 

~ di ~ dX ~ dx dX 
du(t,x{t,X))- 



dx 



-F(t,X) = Vu • F, 
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which is exactly equation (jl.lc|) . Here, and in what follows, we use the conventional 
notations: 

du ■ dF ■ ■ 

(Vu)jj = (VuF)y = (Vu) ik F kj , (U • VF)ij = U k-Q^, 

and summation over repeated indices will always be well understood. In viscoelastic fluids, 
(jl.lcp can also be interpreted as the consistency of the flow maps generated by the velocity 
field u and the deformation gradient F. 

The difference between fluids and solids lies in the fact that, in fluids, such as Navier- 
Stokes equations [26J , the internal energy can be determined solely by the determinant part 
of F (equivalently the density p, and hence, ()l.lcf) can be disregarded); while in elasticity, 
the energy depends on all information of F. 

In the continuum physics, if we assume the material is homogeneous, the conservation 
laws of mass and of momentum become [7J \16\ [31] 

d t p + div(pu) = 0, (2.3) 

and 

d t (pu) + div(pu <g> u) - pAu + VP(p) = div((det F)~ 1 5F T ), (2.4) 
with divu = where 

pdetF = l, (2.5) 

and 

dW 

Here S, pSF T , W(F) denote Piola-Kirchhoff stress, Cauchy stress, and the elastic energy 
of the material, respectively. Recall that the condition f|2.6|) implies that the material is 
hyperelastic [23] . In the case of Hookean (linear) elasticity [151 EH [20] , 

W(F) = i|F| 2 = itr(FF T ), (2.7) 

where the notation "tr" stands for the trace operator of a matrix, and hence, 

S{F) = F. (2.8) 

Combining equations (|2.1|) - (|2.8p together, we obtain system (jl.ip . 
If the viscoelastic system (jl.ip satisfies 

div(p F([) = 0, 

initially at t = with Fo = / + Eq, it was verified in [18] (see Proposition 3.1) that this 
condition will insist in time, that is, 

div(p(t)F(t) T ) = 0, for t > 0. (2.9) 

Another hidden, but important, property of the viscoelastic fluids system ([l.ip is con- 
cerned with the curl of the deformation gradient (see [15, 16J). Formally, the fact that the 
Lagrangian derivatives commute and the definition of the deformation gradient imply 

XkFlJ = dx k dXj = dx 3 dx k = x > Fik > 

which is equivalent to, in the Eulerian coordinates, 

FjfcVjFtfCt, x{t, X)) = F nj V n F ik (t, x{t, X)), 
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that is, 

FjjfcVjFy(i,sc) =F nj V n F ik (t,x), 
which means that, using F = / + E, 

VkEij + EikViEij = VjEik + E n jV n Eik- (2-10) 

According to (|2,10p , it is natural to assume that the initial condition of E in the viscoelastic 
fluids system (jl.ip should satisfy the compatibility condition 

V fe #(0)y + £(0)i fc VjE(0)y = VjE(0) ik + E(0) nj V n E(0) ik . (2.11) 

Finally, if the density p is a constant, (jl.ip becomes its corresponding homogeneous 
(density-independent) incompressible form (see [UlITI^IIElIIBJCEHlEn] and references therein) 

divu = 0, 

<9fU + u- Vu-/jAu + VP = div(FF T ), (2.12) 
<9 t F + u- VF = VuF. 

For more discussions on viscoelastic fluids and related models, see [H [Til HI El 
ES E21 [231 EU El] and the references cited therein. 

3. Main Results 

In this Section, we state our main results. As usual, the global existence is built on the 
local existence and global uniform estimates. 

In this paper, the standard notations for Sobolev spaces W s ' q and Besov spaces Bp q ([3]) 
will be used. Throughout this paper, the real interpolation method ([3]) will be adopted 
and the following interpolation spaces will be needed 

Xf-^ = (L 9 (IR 3 ), W 2 ' q (R 3 )) 1 _ 1 = B^~*\ 

p ' " 

and 



Yp * = {L*$*),W x «<R*)) x _ k = B qP \ 



Now we introduce the following functional spaces to which the solution and initial con- 
ditions of the system (|1.1|) will belong. Given 1 < p,q < oo and T > 0, we set 

Q T = R 3 x (0,T), and 

W p ' q (0, T) := {u G W^ p (0, T- {L q (R 3 )) 3 ) D L p {0, T; {W 2 ' q (R 3 )) 3 ) : divu = 0} 
with the norm 



as well as 



with the norm 



We denote 



l u ll>VP'9(0,T) : ~ II u IIh /1 .p(0,T;L9(R 3 )) + II u IIlp(0,7W 2 .9 

v*« -.= ^''nr^J 3 x (w^iR 3 )) 10 
\\(f,g)\\v*« := H/ll 2(1 _i) + 11/11 li + NIwtm 

y\.p i p 

W(0, T) = W p ' q (0, T) n W 2 ' 2 (0, T), 
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and 

v = v™ n v 2 ' 2 . 

Our first result is the following local existence: 

Theorem 3.1 (Local existence for viscoelastic fluids). Let Tq > be given and (uq, po, Eq) 

£ Vq with p £ [2, oo),q E (3, oo). There exists a positive constant 8q < 1, depending on To, 
such that if 

||(u ,p - l,#o)||v ^ S 0, 
then the initial-value problem (|l.la|) - (|l.lc|) with (|1.2|) as well as the initial-value problem 
(|1.1|) - (|1.2|) have a unique strong solution on M 3 x (0,Tq), satisfying 

(u,p,F) G W(0,T ) x (Ty 1 ' p (0,r ;L 9 (M 3 )nL 2 (M 3 ))nL p (0,ro;Ty 1 '''(M 3 )nVF 1 ' 2 (lR 3 ))) 10 . 

Remark 3.1. The solution constructed in Section 4 later does not require the incompress- 
ible condition (jl.ldp . thus we have the local existence for both the compressible and 
incompressible cases. The solutions in Theorem 13.11 is local in time since 6q = 5q(Tq) 
implies that To is finite for a given 5o 1. 

Remark 3.2. An interesting case is the case q < p. Indeed, by the real interpolation 
method and Theorem 6.4.4 in [3], we have 



and 



?n 1 )a 2 (1- L ) 2 ( 1 - i ) 
W 2[1 p>' q C B qp v =X P ( * J , 



1 1 a 1 - 1 1-- 

W 1 —^ C Bqp p =Z P p . 



Then, if we replace the functional space Vq' 9 in Theorem 13. II by 

V™ := ((Ty 2{1 "p ) ' <? (M 3 )) 3 n(VF 1_ p' 9 (R 3 )) 3 ) x {W^iR 3 )) 10 , 

Theorem 13.11 is still valid. 

The above local existence, with the aid of global estimates and the suitable choice of 
the smallness of the initial data, will result in the following global existence: 

Theorem 3.2 (Global existence for viscoelastic fluids). Assume that 

• p = 2 and q £ (3, 6]; 

• There exists a 5q > 0, such that, for any 5 with < S < So <C 1, the initial data 
satisfies 

||(u ,po - l,^o)||vb < (3-1) 
with compatibility condition (|2.1ip , divuo = 0, and 

div( P0 F r ) = 0; (3.2) 

• In addition, the initial data satisfies 

[ (^p \u \ 2 + lp \E \ 2 + ^—(pl- 1Po + 7 -l))dx<5\ (3.3) 

and 

||Vu || L 2 + ||u • Vu || L 2 + ||Au || L 2 + \\Vp \\ L 2 + HV^olU^ < J 4 . (3.4) 
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Then, there exists a /Uq > depending only on q ( and determined by (|6. 18j) ), such that if 
< fi < Ho, the initial-value problem (jl.ip - ([1.2p has a unique strong solution defined on 
R 3 x (0, oo) with 

(u, p, F) G W(0, T) x (^ 1 ' p (0, T; L 9 (M 3 ) n L 2 {R 3 )) n L p (0, T; W 1,q (R 3 ) n W 1,2 (R 3 ))) 10 , 
/or eac/i T > 0. Furthermore, the solution satisfies 

sup ||(u(i),p(i)-l,£(i))|| yr < V^. (3.5) 
te[o,oo) 

Remark 3.3. Notice that if g > 3, then by Theorem 5.15 in [I], the imbedding Ty 1,l? (IR 3 ) «— > 
C^(IR 3 ) is continuous. Here, the notation C^(M 3 ) means the spaces of bounded, continuous 
functions in R 3 . Hence the condition (|3.ip implies that, if we choose 5 sufficiently small, 
by Sobolev's imbedding theorem, there exists a positive constant Co such that 

Po > C > 0, for a.e. x G M 3 . (3.6) 

Remark 3.4. Under assumption (|3.2|) . the authors in |171 [T8] showed that the property 
will insist in time, that is, for all t > 0, 

div(pF T ) = 0. 

Remark 3.5. If the density p is a constant, for simplicity, p = 1, Theorems 13.11 and 13.21 
become the analogous results of the homogeneous incompressible viscoelastic fluids (|2,12p . 
In other words, following our argument in this paper, we can recover the global existence of 
strong solutions, or even classical solutions, of the homogeneous incompressible viscoelastic 
fluids near its equilibrium. 

An important consequence of Theorem 13.21 is the case as E = and disregarding the 
equation (jl.lcp . In this case, one has the global existence of density-dependent incom- 
pressible Navier-Stokes equations, since the term on the right-hand side of (jl.lbp can be 
incorporated into the pressure. We state the result without proof as follows. 

Corollary 3.1 (Global existence for Navier-Stokes equations). Assume that 

• p = 2 and q G (3, 6]; 

• There exists a 5q > 0, such that, for any 5 with < 5 < 5o <C 1, the initial data 
satisfies 

||(u ,p -l)lk <<^, 

with divuo = 0; 

• In addition, the initial data satisfies 

[ (^Po|u | 2 + ^-(p2:-7Po + 7-l))^<^ 4 - 
7r3 \2 7-1 ) 

Then, the initial-value problem for the density- dependent incompressible Navier-Stokes 
equations has a unique strong solution defined on M 3 x (0, oo) such that 

(u, p) G W(0, T) x (W 1,p (0, T; L q (R 3 ) n L 2 (R 3 )) n L p (0, T; W 1 ' q (R 3 ) n W 1 ' 2 (R 3 ))) , 

for each T > 0. Furthermore, the solution satisfies 

sup \\{u(t),p(t)-l)\\ vr < v 7 ^- 

ie[0,oo) 
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Remark 3.6. The similar result to Corollary 13. II was shown in [§]. It is worthy of noticing 
that there is no assumption on the amplitude of the viscosity \i or condition (|3.4p because 
those two conditions are only useful when dealing with the dissipation of the transforma- 
tion gradient F. Actually, as seen from the argument later in this paper, Corollary 13. II still 
holds if p > 2 and q > 3. 



4. Local Existence 

In this section, we prove the local existence of strong solution in Theorem 13.11 To this 
end, we introduce the following new variables by scaling 

s := vH, y := ux, v(y, s) := -u(x, t), r(y, s) := p(x, t), G(y, s) := E(x, t), 

v 

where v > will be determined later. Then, system (11. ip . with (jl.lcp replaced by (II. 3|) . 
becomes 

rt + div(ry) = 0, (4.1a) 
(rv) t + div (rv ® v) - fiAv + i<- 2 VP = v~ 2 div (r(I + G)(I + G) T ) , (4.1b) 

G t + vVG = VvG + Vv, (4.1c) 
divv = 0. (4.1d) 

From (|2.10p . one has 

VfcGjj + GikViGij = VjGik + G n jV n Gik- (4.2) 

Thus, if we denote by G{ the i-th row of the matrix G (or the i-th component of the vector 
G), then P~2|) becomes 

curl Gi = G ni V n G ik - G lk ViG iy (4.3) 

The proof of local existence of strong solution with small initial data will be carried out 
through three steps by using a fixed point theorem. Instead of working on (jl.ip directly, 
we will work on (|4.ip . We note that (|4.ip is just a scaling version of (jl.ip . It can be seen 
from the argument below that we only need to verify the local existence in W p ' q (0, T), 
< T < To, while initial data belongs to Vq' q . 

4.1. Solvability of the density with a fixed velocity. Let Aj(x,t), j = l,...,n, be 
symmetric m x m matrices in M n x (0, T), f(x,t) and vq(x) be m-dimensional vector 
functions defined in 1" x (0,T) and M n , respectively. 
For the following initial-value problem: 

n 

d t v + A j( x > t ) d i v + B ( x > t) v = f( x > *)> { . 4 n 

i=l ^ ' 

v(x,0) = v (x), 

we have 
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Lemma 4.1. Assume that 

A, e [C(0, T; H s (R n )) n C x (0, T; ff- 1 ^"))] mxm , j = 1, n, 

B eC((0,T),H s -\R n )) mxm , f eC((0,T),H s (M n )) m , v eH s (R n ) m , 
with s > ^ + 1 is an integer. Then there exists a unique solution to (|4.4p . i.e, a function 

v e [C([0,T),H s (R n ))nC 1 ((0,T),H s ' 1 (R n ))] m 
satisfying (|4.4p pointwise. 

Proof. This lemma is a direct consequence of Theorem 2.16 in |26| with Ao(x,t) = /. □ 

To solve the density with respect to the fixed velocity, we have 

Lemma 4.2. Under the same conditions as Theorem \3.1\ there is a unique strictly positive 
function 

r := S{v) e W l ' p (Q, T; L q (R 3 )) n L°°(0, T; W^iM 3 )) 

which satisfies the continuity equation ()4.1ap andr — 1 £ L°°(0, T; L 9 (M 3 )). Moreover, the 
density satisfies 

\N r \\L°°(0,T;Li(WS)) < C(T, |M| W(0 ,T)) ( || Vr 1| L g( R 3) + l) , (4.5) 

and the norm \\S(v) — l\\ w i, q m3\(t) is a continuous function in time. 



Here, and in what follows, C stands for a generic positive constant, and in some case, 
we will specify its dependence on parameters by the notation C(-). 

Proof. For the proof of the first part of this lemma, we refer the reader to Theorem 9.3 in 
[26J, or the first part of the proof for Lemma 14.31 below. The positivity of density follows 



directly from the observations: by writing ([4. lap along characteristics 4rX(t) = v, 

j t r(t, X(t)) = -r(t, X(t))dwv(t, X(t)), X(0) = x, 
and with the help of GronwalFs inequality, 



(inf /? )exp I - / ||divy(t)|| LO o (]R 3)da; ) < r(t 



o 



< (suppo)exp / \\divv(t)\\ Lao r R 3)dx 







Now, we can assume that the continuity equation holds pointwise in the following form: 

dtr + rdivv + v ■ Vr = 0. 
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Taking the gradient in both sides of the above identity, multiplying by \Vr\ q ~ 2 Vr and 
then integrating over M 3 , we get, by Young's inequality 

- — \\Vr\\ q Tafm3 , < f \Vr\ q \divv\dx + f r\ Vrl^ 1 ) Vdivvldx 
qdt u " i9 ( R ) J R3 7 R3 

+ / \Vv\\Vr\ q dx - - f vV\Vr\ q dx 

^ H Vr llM(R3) (IIVuIUoo^rs) + lkllL-»(R3)||Vdivt;|| LlJ(IR 3 ) ) ( 4 - 6 ) 

+ - J divv\Vr\ q dx + ||r|| ioC ( R 3)|| Vdivv || i9( - R 3) 
Q Jm. 3 



< C\\Vr\\l 3 \\v\\ W 2, q(m + ||r|| i0 o (R 3 ) ||Vdivw|| i 
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since q > 3. Then (|4.5|) follows from Gronwall's inequality. 

Finally, noting from P~6j) and (|4"3|) that ^||Vr||^ ?(R3) G L 1 (0,T), and hence 

|«V(r- 1)|IL (M3) €£ l (0,n 

which together with (j4.5j) implies that ||V(r — l)||| g , K3 j(i) is continuous in time, and hence, 
||V(r — l)||is(K3)(t) is continuous in time. Similarly, from the continuity equation, we know 
that 

d t (r-l) = -div((r - l)v) - divv G L p (0, T; L q (R 3 )), 
which, together with the fact r - 1 G L°°(0, T; L q (R 3 )), yields r — 1 G C([0, T]; L q (R 3 )). 
Hence, the quantity \\r — l\\w 1 <i(R 3 )(t) is continuous in time. The proof of Lemma 14.21 is 
complete. □ 

4.2. Solvability of the deformation gradient with a fixed velocity. Due to the 

hyperbolic structure of (|4.1c|) , we can apply Lemma 14.11 again to solve the deformation 
gradient G in terms of the given velocity. For this purpose, we have 

Lemma 4.3. Under the same conditions as Theorem \3.1\ there is a unique function 

G := T(v) G W l ' p {Q, T; L q {M 3 )) Pi L°°(0, T; W 1,9 (R 3 )) 

which satisfies the equation ([4. lc[) . Moreover, the deformation gradient satisfies 

l|VG|| ioc( o !T;L9(R3)) < C(T, \\v\\ w(0)T) ) (l|VG(0)|| L9(R 3) + l) , (4.7) 

and, the norm \\G\\wi,qma){t) "is o, continuous function in time. 

Proof. First, we assume that v G C 1 (0, T; (IR 3 )), G G G^°(1R 3 ). Then, we can rewrite 
(|4.1cp in the component form as 

a t Gj + v ■ VGj = Vv Gj + Vvj , for all 1 < j < 3. 

Applying Lemma [4 . 1 1 successively with Ak(x, t) = Vk(x, t)I for all 1 < k < 3, B(x, t) = Vv, 
and f(x,t) = Vvj, we get a solution 

G G {C^T, fl" i_1 (K 3 )) n {0,T;H l {R 3 ))} , 
which implies, by the Sobolev imbedding theorem, 

G G n^C^O, T; C k (R 3 )) = C x (0, T; G°°(1R 3 )). (4.8) 
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Next, for v S W p,q (0, T), there are two sequences: 

v n eC\0,T;C^(R 3 )), Gq G C^°(IR 3 ), 

such that 

v n -> u in W(0,T), GJ-> G in ^(R 3 ), 

thus w n — > u in C(B(0, a) x (0, T)) for all a > where -B(0, a) denotes the ball with radius 
a and centered at the origin. According to the previous result (|4.8|) . there are a sequence 
of functions {G„}^° =1 C C^O, T; C°°(M 3 )) satisfying 

+ v n • VG n = Vv n G n + Vv n , (4.9) 

with G n (0) = Gft. Multiplying ([49]) by |G n | <7 ~ 2 G n , and integrating over R 3 , using inte- 
gration by parts and Young's inequality, we obtain, 

~ / \G n \«dx 
q at J R 3 

= -- f v n -V\G n \ q dx + [ Vv n \G n \ q - 2 G 2 n dx+ [ Vv n \G n \ q ~ 2 G n dx 

Q JR3 JR3 Jm? 

< — ^-\\G n \\ q Lq (\\Vv n \\ L ^ + ||Vw n || w ) + ||Vu„||m- 

From Gronwall's inequality, one obtains, 

/ \G n \ q dx 
Jm. 3 

/ \G n (0)\ q dx + q \\Vv n \\ Lq exp(- (q+l)(\\Vv n \\ L ^ + \\Vv n \\ Lq )dT)ds) 
Jm. 3 Jo V Jo J J 



x exp ^ (q+ l)(\\S7v n \\L°° + \\Vv n \\ Lq )ds^ 
< ^ \G n {0)\ q dx + q WVvnWitdaJ exp ^jf (g + l)(||Vu n |U<» + ||V^|| W )^) • 
Thus, 

l|G'n||Loo( 0)T;L g( R 3)) < C(T, ||w||lp(0,T;H/ 2 .9(M 3 ))) (II G (0) \\lv(m?) + l) < oo. (4.10) 
Hence, up to a subsequence, we can assume that the sequence {v n } was chosen so that 

G n ^G weak-* in L°°(0, T; L q (R 3 )). 
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Taking the gradient in both sides of (|4.9p . multiplying by |VG n | 9_2 VG n and then 
integrating over K 3 , we get, with the help of Holder's inequality and Young's inequality, 

' "VG n ||* 3 > < / \VG n \«\Vv n \dx+ [ IGnWVGnl^VVVnldx 

V ; JR 3 JR 3 



q dt 



+ j \Vv n \\VG n \ q dx - - I r„V\YG„\''fl.r 



+ f \VVv Tl \\VG Tl \ q - l dx 

< I \VG n \ q \Vv n \dx + I |G n ||VG n | 9 - 1 |VVv n |^ uu] 
Jrz Jr 3 K^- 11 ) 

+ / \Vv n \\VG n \ q dx + - / \X7v n \\X7G n \ q dx 
Jr 3 Q Jr 3 

+ f |VVu n ||VG n | 9_1 dx 
Jr 3 

< C , ||VG n ||^^ R3 ^||i; n || l y2,g(] R 3) + (HGnll^^a^ + 1)11^11 w 2 -« 

< C , ||VG n ||^^ R3 ^||i; n || H /2,g( R 3) + (||G n ||^^ R3 ^ + 1)11^11^/2, 9 ( 

since q > 3. Using Gronwall's inequality and (|4.1U|) . we conclude from (|4.11|) that 
WVG n \\ Loa (o iT . L9 Qn3^ < C(T, ||tVi||yv(o,T)) (H v Gn(°)llL9(R 3 ) + 1) > 

and hence, 

ll^ 7 G'|| L oo( 0jT;L « 7 ( K 3)) < liminf ||VG n || L oo( 0jT;i ,( R 3)) 

<C7(T 1 ||t;|| w(0)r) )(||VG(0)|U, (R 8 ) + l). 
Therefore, 

l|G ! n||L°o(0,T;W/ 1 .9(R3)) < C(T, \\v ||lp(0,T;W 2 >«(R 3 )) > II ^(°) II W^R 3 )) < 00 • 

Furthermore, since q > 3, we deduce G S L°°(Qt) and 

\\G\\l°°(Q t ) < C(T, ||-y||iP(o,T;W 2 .9(R 3 ))) 11^(0)11^/1, 9 ( K 3)) < OO. 

Passing to the limit as n — > oo in (|4.9p . we show that (|4.1cj) holds at least in the sense of 
distributions. Therefore, 

d t G £ L p (0,T;L q (R 3 )), 



(4.12) 



then G G W^ p (0,T; L q (R 3 )), and hence G 6 C([0,T];L q ( 

Finally, to show that the quantity ||Cr||iyi,?(R3)(i) is continuous in time, it suffices to 
show that || VG||x,,m3) is continuous in time. Indeed, from (|4.11|) . we know that 

f t \\VG\\l qim (t)eLno,T), 

which, with (|4.12j) . implies that ||VG|| i9 ( R 3) G C([0, T]). The proof of Lemma 14.31 is 
complete. □ 
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4.3. Local existence via the fixed-point theorem. In order to solve locally system 
(|4.ip . we need to use the following fixed point theorem (cf. 1.4.11.6 in |26j): 

Theorem 4.1 (Tikhonov Theorem). Let M be a nonempty bounded closed convex subset 
of a separable reflexive Banach space X and let F : M i— > M be a weakly continuous 
mapping (i.e. if x n G M,x n — ► x weakly in X, then F(x n ) — ► F(x) weakly in X as well). 
Then F has at least one fixed point in M . 

Now, let us consider the following operator 

Lu j:=^- /jAuj, uj G W p ' q (0, T). 

One has the following theorem by the maximal regularity of parabolic equations; see 
Theorem 9.2 in [26], or equivalently Theorem 4.10.7 and Remark 4.10.9 in [2] (page 188). 

Theorem 4.2. Given 1 < p < oo, ui G V^ q and f G L p (0, T; L 9 (M 3 ) 3 ), the Cauchy 
problem 

Lio = f, i€(0,T); w(0)=w o , 
has a unique solution uj := L~ x (u)q, f) G W p ' q (0,T), and 



|wp>9(0,T) < C ^[|/[|iP(Q,T;L9(R3)) + 1 1 ^0 1 1 V P ) : 



where C is independent of ujq, f and T. Moreover, there exists a positive constant cq 
independent of f and T such that 



IMIw^(0,T) > CO SUp \\UJ{t)\\ V P,i. 

te(o,T) 

Notice that Theorem 14.21 implies that the operator L is invertible. Thus we define the 
operator H(y) : W p ' q (0,T) h+ W p ' q (0,T) by 

U{y) := L- 1 (v , d t {(l - S{v))v) - div{S(v)v ® v) + u- 2 V{P{l) - P(S(v))) 

+ v- 2 div(S(v)(I + T(v))(I + T(w)) T )) . (4.13) 

Then, solving system (14. ip is equivalent to solving 

v = H(v). (4.14) 

To solve (I4.14p . we define 

B R (0) := G W M (0,T) : |b|| WP , 9(0 , T) < i?}. 
Then, we prove first the following claim: 
Lemma 4.4. There are v, T > 0, and < R < 1 such that 

H(B R (0)) C B R (0). 

Proof. Let T>0, 0<i?<l and v G Br(0). Since <S(u) solves (|4.1aj) . we can rewrite 
operator Ti as 

U{v) = L- 1 fv , (1 - S{v))d t v - S{v)v ■ Vv + ^ 2 V(P(1) - P(5(«))) 

+ ^- 2 div(S(t>)(/ + T(v))(I + T{v)) T )) . (4.15) 
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Thus, it suffices to prove that the terms in the above expression are small in the norm of 
LP(0,T; (L<?(K 3 )) 3 ). 

First of all, we begin to deal with the first term by letting r := S(v) — 1. Thus, r 
satisfies the equations 

j dtf + div(fv) = 0, 
} r(x, 0) = ro — 1. 
Repeating the argument in Section 14.21 again, we obtain 

II^IU^CQt) - C|I^IIl oo (0,T;VK 1 '«(K 3 )) 

< C \\ r - l||wl.?(M3)C( r , IMIw(0,T)) 

< \\r -l\\ w i, q{M3) C(T,R)<C(T)R, 
where, by the formula of change of variables, we deduce that 

3 

ll r - M\li(R3) < VI ||/9 - 1||l<J(R3) < R, 



and 

3-q 

l|Vr ||z,?(K3) < v " ||V/9o||l8(r3) < R, 
if \\po — 1\\l<!(r 3 ) is small enough and v > 1 is large enough. Hence, due to the assumption 
v 6 -Br(O), we obtain 

||(1 - S(v))d t v\\ LP ^ T , Lqm) < C(T)R 2 . (4.16) 
Secondly, by the Sobolev imbedding, 

\vVv\ q dx < \\v\\ q Lq(M 3)\\Vv \\ L<x m3) < C|l u llLq(K3)ll^ll^/2, 9 



and thus, since W 1,p (0, T; L q (R 3 )) ^ C([0, T];L q (R 3 )), we deduce 

T / r \ f rT 

rY)'\' l fl:r ) ds < C J \\v\\ P Lq(R3) IMI^/2, 9(m 3 

< C\\v\\ P r x(Vl _T:r 1 i(M'i))\\ v \ lP 



(4.18) 



lL°°(0,T;L9(R 3 )) II u "LP(0,T;W 2 >1 

Therefore, we get 

\\S(v)(v ■ V)v\\ LP{0tT . Lqm) < CR 2 . (4.17) 
Thirdly, for the term VP(5(u)), we can estimate it as follows 

< C(T) sup {P'( V ) : C{T)~ l < V < C(T)} (||Vr |U 9( R 3) + 1) 
Fourthly, for the term div(S(v)(I + T{v)){I + T(v)) T , we have 
|div(S(t;)(/ + T(v))(I + T(v)) T )\ < \VS(v)\\I + T(v)\ 2 + 2S(v)\ VT(v)\ \I + T(y)\, 
and hence, 

||div(5(«)(/ + T(v))(I + T(v)) T )\\ LP{0 ^ Lq) 

< 111^^(^)11^ -h -^(^)| 2 ||^co,Gr;i.) 2||^(^)| + ^(^)| ||^ ( o,^;^) (4-19) 

< C{T)M, 
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with 

M = max{||G ||w-i,« + 1, ll^o ||z,°°(R 3 ) + h ll r o||iyi><? + 1, IMIl°°(r3) + l} 3 < oo. 

Combining together (|4.16|) , (|4.17|) , (|4.18|) , (|4.19|) , using the Theorem 14.21 and assuming 
parameter v sufficiently large and R < 1 sufficiently small, we get 

\\H{v)\\ w ^ T) <C(T)(R 2 + v- 2 )<R. 

The proof of Lemma 14,41 is complete. □ 

Thus, it is only left to show the following: 

Lemma 4.5. The operator TL is weakly continuous from W p ' q (0, T) into itself. 

Proof. Assume that v n — > v weakly in W p ' q (0, T), and set r n := S(v n ), G n := T(v n ), then 
{r n }n=i and {G n }%=i are uniformly bounded in L°°(0, T; W 1 ' 9 (R 3 )) n W 1,p (0, T; L q (R 3 )) 
by Lemmas 14. H and 14. 31 Hence, up to a subsequence, we can assume that r n — > r and G n — > 
G weakly* in 7°°(0, T; W 1 ' q (R 3 )) n W ljp (0, T\ L q (M. 3 )) and then strongly in C((0,T) x 
73(0, a)) for all a > 0. And at least the same convergence holds for v n . Thus, (|4.1a)) and 
(|4.1cp follow easily from above convergence. 

Since r n -» r weakly* in 7°°(0, T; W 1 ' q (M 3 )) n W 1,p (0, T; 7 9 (M 3 )), we can assume that 
P'(S(v n ))VS(v n ) -» P'(<S(u))V<S(u) weakly in 7 P (0, T; 7 9 (IR 3 )) and hence, 

7" 1 (0, VP(5(v„))) -» L" 1 (0, VP{S(v))) weakly in W(0, T), 

since the strong continuity of 7 _1 from 7 P (0, T; 7 9 (R 3 )) into W(0, T) and the linearity of 
the operator L imply also the weak continuity in these spaces. 

Similarly, since d t v n -> d t v weakly in L p (0, T; L 9 (M 3 )) and r n -» r in C((0, T) x B(0, a)) 
for all a > 0, we have (r e — r n )dtv n — ► (r e — r)9fW weakly in 7 P (0, T; 7 9 (M 3 )) and conse- 
quently 

L- 1 (0, (r e - r„)d t <) 7" 1 (0, (r e - r)d t v) weakly in W(0, T). 

Since Vv n weakly in W^O, T; PF _1 ' 9 (R 3 )) n 7 P (0, T; VF 1 ' 9 ^ 3 )) which is com- 

pactly imbedded in to C([0, T]; L q (B(0, a))) for all a > 0, we can assume that w n — > i> 
strongly in 7°°(0, T; L q (B(0, a))) for all a > 0, and then 

<S(«n)( w n • V)f„, -> <S(v)(u • V)t> 

weakly in 7 P (0, T; 7 9 (M 3 )). Hence 

2T 1 (0, 5(u n )(u n • V)u„) -» 7 _1 (0, • V)«) weakly in W(0, T). 

Finally, due to the facts that r n -> r and G n -> G weakly* in 7°°(0, T; W 1 ' 9 ^ 3 )) n 
W^O.TjL^fR 3 )) and strongly in C((0,T) x B(0, a)) for all a > 0, we deduce that 

div(5(v n )(J + TK))(/ + T(v n )) T ) - div(5(«)(J + T(i;))(7 + 7») T ) 
weakly in L p (0, T; L 9 (M 3 )). Therefore, 

L- 1 (o, div(5(v n )(J + T(v n ))(I + TK)) T )) -» L- 1 (o, div(5(v)(7 + T(«))(J + T(t;)) T )) 

weakly in W(0,T). 

Thus, we can conclude that 

H(v n )^H(v) weakly in W(0,T), 

due to the weak continuity of map 7 _1 (w, 0). The proof of Lemma 14.51 is complete. □ 
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Therefore, by Theorem 14.11 there exists at least a fixed point 

v = H(v) G B R (0) C W(0, T), (4.20) 

of mapping H. The fixed point v provides a local in time solution (p, u, F) of system (jl.ip 
near its equilibrium through the scaling with v sufficiently large. 

The proof of the local existence in Theorem 13.11 is complete. The uniqueness will be 
proved in the next section. 



5. Uniqueness 

In this section, we prove the uniqueness of the local solution found in the previous 
section. Notice that, the argument in Section U] yields that 

d t v G L 2 (0,T;L 2 (M 3 ), Vr G L 2 (0, T; L 2 (M 3 )), VG G L 2 (0, T; L 2 (R 3 )). 

Hence, using the interpolation, we deduce that 

dtv G L P0 (0,T;L 3 (M 3 ), Vr G L Po {0, T; L 3 (R 3 )), VG G L Po (0, T; L 3 (M 3 )), 

where 

i _ e l-e i _ e i-e 

Po 2 p ' 3 2 g 
for some # G [0, 1]. Now, assume that v\, v-i satisfying (|4.20p for some T > 0. Let 

r := S(vi) -S(v 2 ), v:=vi-v 2 , G := T(vi) - T(v 2 ), 

with a little abuse of notations (however, there should be no confusion in the rest of this 
section). Then, we have 

dtr + v i • Vr + v ■ VS(v2) + rdivf i + S(v2)dwv = 0, 

r(0) = 0. " (5 - 1} 



Multiplying (|5.ip by r, and integrating over M 3 , we get 

__ rll r llf2 / \r\ 2 d\vv\dx+ I ( vV S (v2)r + \r\ 2 divv\ + rS (v2)divv) dx = 0, 

2 at 2 j R 3 j R 3 

which yields 

j t \\r\\ 2 L 2 m < ||div7; 1 || L ^||r|| 2 2 +e||V«|| 2 2 +C7(e)||VcS(z; 2 )r|| 2 f 

+ e\\Vv\\l 2m + C(e)\\S(v2)\\ 2 LO o\\r\\ 2 2 
< ||div^ 1 || ioo ||r||| 2 +£||V U || 2 :2 +C( e )||V5(^)|| 2 :3 ||r||| 2 ( 5 - 2 ) 

+ 4^v\\h m + C(e)\\S(v2)\\ 2 L o a \\r\\ 2 
< m (e)\\r\\ 2 L2 +2e\\Vv\\ 2 L2( 



(5.3) 



where e > 0, ^(e) = ||divui||Loo + C(e) (||V5(^)||» 3 + ||5(«2)|||oo) . 
Similarly, from (|4.1cj) . we obtain 

J d t G + v 1 -VG + vVG2 = Vv 1 G + VvG 2 + Vv, 
1 G(0) = 0. 
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Multiplying (|5.3p by G, and integrating over R 3 , we get 



ld_ 
2~dt 



G \\\z ~lf \G\ 2 div Vl dx + f v ■ VT{v 2 ) : Gdx 

2 J R 3 J R 3 

/ \G\ 2 Vv 1 dx + [ VvT(v 2 ) : Gdx + / : Gdx, 



which yields 
d 



\\G\\i HR3) < lldiv^lliocHGII^ +e\\Vv\\b + C{e)\\VT{v 2 )G\\\ 

+ e||V^||| 2(M3) + C(e)\\T(v 2 )f LX ||G||| + e||V<;||| 2 + G( £ )||G||| 2 
< lldiv^ll^HGII^+ellVHll.+C^llVTMIllallGlll, (5-4) 

+ 4Vvf L2m + C(e)\\T(v 2 )\\ 2 Lx \\G\\l + e\\Vv\\ 2 L2 + C(e)\\G\\ 2 L2 
<V2(s)\\G\\ 2 L2 +3e\\Vv\\ 2 L2( 



where 772(e) = ||divwi||^ + C(e) (\\VT(v 2 )\\ 2 L3 + II^MUloc + l) . 
For each Vj, j = 1, 2, we deduce from (|4.1bj) that 

S(vj)dtVj — fiAvj 

= -SfaXv, ■ V)vj - VP(S(vj)) + div(5(^)(7 + T( Vj ))(I + T( Vj )) T ), 
k Vj(0)=v . 

Subtracting these equations, we obtain, 

S(vi)d t vi - S(v 2 )d t v 2 - fj,Av 

= -<S(^)(^ • V)«i + SfaX^ • V)v 2 - VP(S(«i)) + VP(S(t> 2 )) (5.5) 
+ div(5(t;i)(/ + T(«i))(/ + T^x)) 7 ) - div(S(u 2 )(J + T(v 2 ))(J + T(v 2 )) T ). 

Since 

- «S(ui)(ui • V)vx +S(v 2 )(v 2 ■ V)v 2 

= -S{ Vl ){v ■ V) Vl - (S(vi) - S{v 2 ))(v 2 • V)«i - S(v 2 )(v 2 ■ V)v, 

and 

5(vi)(J + T( Vl ))(I + T( Vl )) T - S(v 2 )(I + T(v 2 ))(I + T(v 2 )) T 
= S( Vl )G(I + T( Vl )) T + r(I + T(v 2 ))(I + T( Vl )) T + S(v 2 )(I + T(v 2 ))G T , 

we can rewrite (15.51) as 



S(vi)dtv — fiAv 

= -rd t v 2 - S( Vl )(v ■ V)«i - 5(tj 2 • V)«i - 5(tj 2 )(tj 2 • V)v - VP(S(«i)) + VP(S{v 2 )) 
+ dw(S( Vl )G(I + T( Vl )) T + r(I + T(v 2 ))(I + T( Vl )) T + S(*, 2 )(J + T( V2 ))G T ). 

(5.6) 
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Multiplying (|5.6p by v, using the continuity equation (jl.laj) and integrating over R 3 , we 
deduce that 



™ / S( Vl )\v\ 2 dx+ [ M |V 



v\ 2 dx 



-S{v\)(v\ ■ V)v ■ v — rdtV2V — S(v\)(v ■ V)v±v — S(v 2 ■ V)i>ii> 

- S(v 2 )(v 2 ■ V)vv - VP(S(vi))v + VP(S{v 2 ))v 

- + T( Vl )) T + r(J + T{v 2 ))(I + T( Vl )) T + S(v 2 )(I + T(v 2 ))G T )Vvdx 

e||v^ni 2 + c( e )|| t s(^ 1 )lli-ll^illi-ll^lli 2 + ^liv^lli^ + ^(^)ll£ > ^2iii 3 |kiii 2 

+ ||S(m)|Uoo||V«i||L<»||u||£a + 2 \\v 2 \\ L °o\\Vv 1 \\ L °o(\\rf L 2 + \\v\\ 2 L2 ) 
+ e||V«||i 2 + C( E )||5(^)||i«||^||ioo||t;||i 2 +e||Vt;||i a 
+ C(e)(mp{P'( V ) : C(T)- 1 < v < C(T)}) 2 ||r|| 2 2 + e||V<;|| 2 2 
+ G( e )(||5(« 1 )||i 00 (l + ||r(« 1 )||io.)||G||i a 

+ 115(^)111^(1 + \\T(v 2 )\\ 2 LOO )\\G\\ 2 L2 + ||r||| 2 (l + ||T(^ 1 )||| 00 )(1 + \\T (v 2 )\\ 2 Laa )) 



(5.7) 



<5e||V ? ;||i 2 +r ? 3(e)(||r||i 2 + || W ||i 2 + ||G||i 2 ) 
with 

m {e) = C&WS^Wl^WvtWlo. + C( £ )||^ 2 ||| 3 + ||5(tn)|| L -||V«i|| L oc 
+ 2||^|| L «||Vt; 1 |Uco + C(e)||5(«2)||ioo||^||i=o 
+ C(e)(snp{P'(r,) : C(T)- 1 < r? < C(T)}) 2 

+ C(e)(\\S(vi)\\l<~(l + \\nvi)\\l°o) + l|S(t»)||L»(l + r(«4)||ioo) 

+ (i + ||r(« 1 )||i.)(i + ||T( U2 )||i«»)). 

Summing up (|5.2p . (|5.4p . and (|5.7| . by taking e = we obtain 
4/ (5(« 1 )|r;| 2 + |r| 2 + |G| 2 )dx + / u / |Vtfdz 

< 2(i»( e ) + r/ 2 ( £ ) + + \\r\\h + \\G\\b) (5-8) 

<2»j(£,i) / (5(^i)|t;| 2 + |r| 2 + |G| 2 )^, 
Jr 3 

with 

+ 772(e) +??i(e) 



?7(e,t) 



min{ min^gRa 5 (t>i ) (x, t) , 1} 

It is a routine matter to establish the integrability with respect to t of the function rj(e, t) 
on the interval (0, T). This is a consequence of the regularity of V\,v 2 € W(0,T) and the 
estimates in Lemmas 14. 21 and 14. 31 for S(vi), T(vi) with £ = 1, 2. Therefore, (|5.8p . combining 
with Gronwall's inequality, implies 

/ (SMM 2 + |r| 2 + |G| 2 ) dx = 0, for all t £ (0,T), (5.9) 
it 3 
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and consequently 

v = 0, r = 0, G = 0. 
Thus, the uniqueness in Theorem 13,11 is established. 



6. Global A Priori Estimates 

Up to now, we prove that for any given To, we can find a unique solution to the scaling 
system ()4. 1|) . That is, we have proved the local existence of solution to the viscoelastic 
fluid system (jl.ip and its uniqueness. In order to establish the global existence for the 
unique solution we constructed in the previous sections, we need to obtain some uniform 
a priori estimates which are independent of the time T. To simplify the presentation, we 
will focus on the case v = 1, that is, system (|l.ip . 

We introduce the new variable: 

a := V In p. 

Then, we have 



Lemma 6.1. Function a satisfies 

d t a + V(u • a) = 

in the sense of distributions. Moreover, the norm \\a 



L'l 



Proof. We follow the argument in [26] (Section 9.8) by denoting a, 
standard mollifier in the spatial variables. Then, we have 



(6.1) 

is continuous in time. 

S E a, where S e is the 



d t o- £ + V(u • a £ 



with 



V(u • a e ) - S e V(u ■ a) 
= (u • Va e - S e (u ■ V<t)) + (a £ Vu - S E {a ■ Vu)) 

=: nl + nl. 

Since a G L°°(0,T;L<? 
[26j (cf. Lemma 2.3 in 

|| (a E - o-)Vu|| L i (0iT;L9 



and u 6 L p (0, T; Vy 1,oc (]R 3 )), we deduce from Lemma 6.7 in 
that TZl — » as e — ► 0. Moreover, 



< Ik 



E "iA (0,T;Li 



IVul 



LP(0,T;L° 



0. 



and S E {c ■ Vu) -> cj • Vu in L^O, T; L 9 (M 3 )) since a • Vu G L p (0, T; L 9 (IR 3 )). Thus, we 
have -> in LP(0, T; L q (R 3 )). Then, taking the limit as e -> in ([E2|), we get flEE]) 



Multiplying (|6.ip by |cr| 9 2 cr, and integrating over 1R 3 , we get 



d II 119 

dt M ^ 



'j u k ~j cr k\ ~ 



ig-2 



1 



divulal 9 )^ 



1 



< || Vull^oo \\a\\ Tq H — ||divu||io 
q 

<c\\ U-||^2,qr||o'||^g. 

Dividing the above inequality by HerHf^ 1 , we obtain 



d „ 

dt MLq 



< C\\\l\\ W 2,q \\0-\\l1- 
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Since a G L°°(0, T; L q (R 3 )), |||o-|| L « G LP(0,T). Thus, ||<r|| L , G C(0,T). The proof of 
Lemma |6. II is complete. □ 

For a given i? = 5q <C 1 as in Section^ if the initial data satisfies ||u(0), po— 1, J5(0)[|vb < 
5 2 with < 5 -C min{^, ^o}, let T(i?) be the maximal time T such that there is a solution 
of the equation u = TC(u) in i?/j(0). By virtue of Lemma [4.2l Lemma [4.31 and Lemma 16. 11 
we know that ||«5(u) — l||wi>9(R3), ||c||l<j and ||T(u)||^i,, are continuous in the interval 
[0, T(R)). On the other hand, under the assumptions on initial data and Remark 13.31 we 
know, if 5 is sufficiently small, then 

lk(o)llL*(R3) < ^l|Vp(o)|| L , (R3) < si « i. 

Hence, there exists a maximum positive number T\ such that 

max{||cS(u) -l|| w i,,(t),||ff|| £ ,(t),||T(u)|| w i,,(t)} < v^R < 1 for all t € [0,Ti]. (6.3) 

Now, we denote T = min{T(i?), Ti}. Without loss of generality, we assume that T < oo. 
Since q > 3, we have 

< c\\p - liityi,,^) < cVr < -, 

if R is sufficiently small. Hence, one obtains 

1 3 
- < p < -. 

2 ~ H ~ 2 

On the other hand, for any given t G (0, T), we can write 
l|u(t)||^ = ||u(0)||S,+jf*^||uW||i,(fa 

= Klli. + ? f(Mt)\W 



Q Jo 



' f \u(s)\ q - 2 u{s)d s u{t)dx] dt 
Jm.3 J 

< ll u o|lz, 9 + - / W^^f^Wds^hi ds 
Q Jo 

<S 2 p + p( f \\u\\ p Lq ds) V ( f \\d s uf Lq ds 



Q 

q 



and consequently, 



HL°°(o,t;L«) + < CR, t G (0, T). (6.4) 



Similarly, we have, for all t G [0, T], 



l u llL°°(0,t;L 2 ) < CR. 
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6.1. Dissipation of the deformation gradient. The main difficulty of the proof of 
Theorem 13.21 lies in the lack of estimates on the dissipation of the deformation gradient. 
This is partly because of the transport structure of equation (jl.lch . It is worthy of pointing 
out that it is extremely difficult to directly deduce the dissipation of the deformation 
gradient. Fortunately, for the viscoelastic fluids system (jl.lj) . as we can see in [Bl [TU 
US [HI CEHl [191 [20], some sort of combinations between the gradient of the velocity and 
the deformation gradient indeed induce some dissipation. To make this statement more 
precise, we rewritten the momentum equation (jl.lb|) as, using (jl.lap 

d t u - (j,Au - divE = -p(u ■ V)u - VP(p) + div(p(I + E) T ) 

+ div((p - l)E) + dw(pEE T ) + (1 - p)d t u, (6 ' 5) 

and prove the following estimate: 
Lemma 6.2. 



Jr 3 

where £ is the fundamental solution of the Laplacian —A in M 3 . Then, (|6.5p becomes 




(6.6) 



Proof. Now we introduce the function Z\{x,t) as 




(6.7) 




(6.8) 



where, with the help of Remark 13.41 

T x = -p(u ■ V)u - VP(p) + div((p - l)E) + div(pEE T ) + (1 - p)d t u. 
Also, from (II. 3D . we have 




(6.9) 



From (16. 8h and (16. 9h . we deduce, denoting Z = u — j^Z\, 

d t Z - pAZ = F := T x - T % 



(6.10) 



where 



?2 



-U+-S* div(Vu£; - (u • V)E). 
fi p 



Equation (|6.10p with Theorem 14.21 implies that 




(6.11) 
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Next, we estimate \\^i\\Lp(o,T;Li(R 3 ))^ i = 1, 2, term by term. Indeed, for T\i using (jlT 
we have 

ll^ r i||LP(0,T;L9(R 3 )) — \\p\\l°°(Q t ) ll U llz°°(0,T;£!(R 3 )) II Vu|| i p(o i Tji°°(R 3 )) 

+ a||V^||xp( 0i T;L8(n)) + O:\W\\lp(0,T;Li{Q,))\\E\\l°°(Q t ) 
+ IIP - 1 II L°°(Q T ) II ^-^11 LP (0,T;Li 



+ H (T llLP(0,T;L'J(R3))||i?||Lo O (Q T ) + || f|| Lp(0,T;L9(R3)) ||-^|||oo (q t ) 
+ 1 1 E 1 1 L°° (Q T ) 1 1 V.E 1 1 LP (0,T;L1 (R 3 )) 

+ \\P~ M\L°°(Q T )\\dtv\\LP(0,T;Li(R 3 )) (6.12) 
< 2^||u|| L p( 0jT;H /2,9( R 3)) + a||V£|| L p( 0jT;i ,j( R 3)) 

Ilp(0,T;L9(R 3 ))I|-E'IIl° o (Q t ) + V /r ^ll cr llLP(0,T;i8( 



+ Vr\\ve\ 



+ a||(j| 

+ R\\ (T \\LP(0,T;L<1QSL 3 )) ^ V ||LP(0,T;£9 

3 / 

< i?2 + a||V-E||£p( 0) T;.W(R 3 )) + V •K||f||iP(0,T;i* 
+ V / ^||V£ , || L p( 0)T;i q( R 3)). 

Here, a = sup {xP'(x) : i < x < |} and in the first inequality, we used the identity 

Vp = -pdivE T - VpE T 

due to Remark 13. 4[ And, for T2, noting that |V£| < Cpjpj an d from integrating by parts, 
we have 

\H < -|u| + --^ * (VuE - (u • V)E), 
p p \x\ z 

with 

\\VuE - (u-V)E\\ 39 < ||Vu|| L p( O T .i3( R 3))||£;|| i oo(o T .i<i 

LP(0,T;L9T^(R 3 )) v v " v 

+ II u IIlp(0,T;L 3 (R 3 )) II V£?||ioo( 0) T;i« 
< i?f . 

Hence, one can estimate, by L p — L 9 estimate of Riesz potential, 



ll^ll 



LP(0,T;Li 



1„ 
< - U 



LP(0,T;H 



c 
+ — 



1 



* (Vu£ - (u • V)E) 



■r 



LP(0,T;L1 



<-R + -\\VuE-(u-V)E\\ 3q 

p p LP{0,T;L^ 

< -(r + crI). 



Therefore, from (|6.12p and (|6.13p . we obtain 



LP(0,T;U 



3 1 

< i?2 + -i? + CK||V-E?||j>fO T-Li 
p ' ' 



+ Vr\m 



LP(0,T;L1(R 3 )) 



+ Vr\\ve\ 



(6.13) 



(6.14) 



LP(0,T;Li 
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Inequalities (|6.1ip and (|6.14p imply that 

\\Z\\lp(0,T;W 2 -i 



< C(p, q) (2R + -R + a\\VE\\ L p^ T . Lq(R3)) + "/R\\cr\\ L p^ T . Lq (^3^ (g _ 15 ) 

+ V^R||V£'||tj>(0 1 T;i9 

Hence, we have, from fl6.Tf) 

HdivEHiP^TjM^ 3 )) 

< M {\\ Z \\lp(0,T;W 2 -i{R 3 )) + II u IIlp(0,T;W /2 .'?(R3))) 

< C(p, q)n(3R + -R + a||V-B|| iP(0)X . L 9( R 3)) + v^IHI^o^.^s)) ( 6 ' 16 ) 



+ v A R||V£ , || L p( 0iT;i q( R 3)) 
On the other hand, from the identity (|4.3p . we deduce that 

1 1 Curl #i|Ui>(0,T;.W(R 3 )) < 2||i?|| L oo(Q T )||Vi?|| L p( 0jT;i9 ( 

< C , ||£ , || L oo( 0iT . H /l, 9 ( R 3))||VS|| L p( 0)T;i9 ( IR 3)) (6.17) 

< CV r R||V-E'||ip(o ) T;M(M 3 ))- 

Combining together (|6.16p and (|6.17p . we obtain 

II v - s IIlp(o,T;L9(r3)) < C(p, q)fi(3R + -R + a|| Vi?|| LP ( 0jT;L q( R 3)) + v^RIHIw ( ,t ; L8 

v n 

+ V^R||V£ , || iP ( 0jT;i q 

and hence, by choosing ^/R <C ^ and the assumption 

C(p,q)pta< 1, (6.18) 

one obtains (|6,6p . The proof of Lemma 16.21 is complete. □ 

Remark 6.1. Notice that, in view of the above argument, estimate (|6,6p is actually valid 
for all t G [0,T], that is, for all t £ [0,T], 



l|V£ , || L p(o, i ;L9(R3)) < C(p,q,n) [R + V R\\a\\ L p( 0AL . 

6.2. Dissipation of the gradient of the density. To make Theorem 13.21 valid, we need 
further the uniform estimate on the dissipation of the gradient of the density. 

Lemma 6.3. For any t £ (0, T), 

lkllLP(o,t;Ls(lR 3 )) ^ C(p,q,fi)R. (6.19) 
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Proof. Multiplying ([Lib]) by o~\o~\ q 2 and integrating over R 3 , we obtain 



a\\ q + / P P\p)\o-\ q dx 



p + A d 

~^Tdi ll0tiLq 

= p Au • a\o-\ q ~ 2 dx — / pdtu ■ a\a\ q ~ 2 dx 
Jm.3 Jm.3 

- [ p(u-V)u- a\a\ q - 2 dx - (p + A) / V(u • a) ■ o-\a\ q ~ 2 

+ [ div(p(I + E)(I + E) T )-a\a\ q - 2 dx. 
We estimate the right-hand side of (|6,20p term by term, 

Au • a\a\ q ~ 2 dx < [|Au||£« [|c[||,9 ; 



< II^uIImIIo-III/; 



(6.20) 



dx 



/ pdtM ■ cr\a\ q 2 dx 
pu • Vu • a\a\ q ~ 2 dx 



< ||U • VU-Hig H^llig 1 

< IIuIImIIuIIw^IMII? 1 ; 

/ V(u • a) ■ a\a\ q ~ 2 dx = / dj\i} t akO'j\a\ q ~ 2 dx+ I u^djd^n p)d ■ In p\a\ q ~ 2 dx 

Jm.3 J K 3 ' J R 3 

I djUkakCj \a\ q ~ 2 dx H — / Ukdk\cr\ 2 \a\ q ~ 2 dx 

JR3 2 J R 3 

/ dj\ikO-kO-j\cr\ 9 ~ 2 dx + / Ufe9fc|o"| \a\ q ~ 1 dx 
/ djUkakCTj\a\ q ~ 2 dx H — f Ukdk\o~\ q dx 

Jr3 q Jr3 

/ djVLkO~kO-j\o-\ q ~ 2 dx I \a\ q divudx 

JR3 Q Jr3 

< C||Vu||z,oo||cr||| a < C||u|| W 2,g \\0-\\ q Lq , 

and, due to (|2,9p . we can rewrite 



div(p(7 + E){I + E) r ) \ . = dA£h±^>h±m 

* 1 CJJb A 



{ei + Ei) d{p(e i + Ej) ) + ^ + Ej) d^±E, 



dxj 



p(ej + Ej 



dEj 
dxf 



then one has 



div(p(J + E)(I + E) T ) • a\a\ q - 2 dx 
< \\VE\\ Lq \\I + E\\ L ~ Halli; 1 < 2\\VE\\ Lq(R3) \\a\\% x . 



dxj 



dE' 

p(ej + Ej)-^-ai\a\ q ^ 2 dx 
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On the other hand, we have 



pP'(p) = P'(l) + (p - 1) / {P'( V (p - 1) + 1) + {flip - 1) + l)P"(v(P " 1) + 1)) d V , 

Jo 



for p close to 1 and consequently 



\\pP'(p) - P'(l)\\ L ~ <\\p- sup ||/(z)| : i < x < i\ 

< C>/Ssup||/(x)| : \ < x < ^\ < CVR, 



where f{x) = P'{x) + xP"(x). Thus, from (16.20p . we obtain 



q dt 

< ||cr||^ 1 ^||Au||L< ! + ||<9 t u|| L 9 + ll u llH/2, q ||u|| L9 

+ C||u|| w a,,||o-|| w + C||V£|| L9(R 3) + ^/R\\a\\ Lq 

< Clklll^f ||Au|| L9 + ||9iu|| L9 + Hullv^HullLa 



+ llullw^lklk 9 + I|V-K||m( R 3) + VR\\a\\ Lq j, 
and hence, by assuming that R <C 1, one obtains 



^ + A d 1 

— «wi, + ^(Dikiii, (621) 

< C 'll cr III? 1 (ll^ u IU« + II^ u IIl" + II u IIw 2 ''j||u||l'J + ||u|| w a,<,||cr||£,, + \\VE\\ Lq ). 



Multiplying (I6T2TT) by || we obtain 

fi -\- X d .. ..p 1 . . p 

< C|kllL9 1 (ll Au llL9 + ll^u|| L9 + ||u|| W 2,,||u|| L9 + ||u|| W 2, 9 ||cr||i g + \\VE\\ L q). 
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Integrating the above inequality over the interval (0, t), one obtains, by using (|6.6 
H + A 1 rt 



P 



\a(t)\f Lq + \p'il)J Q \WWlods 



P-l / I 
p If ft \ p 



+ (HullL^CCtiL^M 3 )) + ll cr ll J L oo (0,i;Z, < ?(R 3 )) + 1 ) / ll U lliy2, 9 ds 
+ \\^ E \\LP(0,t;L"(R' A ))) 

<^^\\a(0)\\ p Lq+ C(p,q,fi)VR\\ 

<T Hlp(0,T;L9(R3)) 

+ C(p,q,n)R[ \\a\\ P Lq ds) (l + llo'llLocCO.fcLafK 3 )) + ll u ll 







lL°°(0,t;L9(R 3 )) "T ll u llL°°(0,t;L9 

and hence, by letting R be so small such that C(p,q, //)vf2 < |, one obtains 
/U + A 1 '* 



\a{t)\\l q + \p'(l) \\a\\l q ds 



p-i 

< ^-^||<T(0)||i 8 + C(p,q,n)R ( f \\a\\ P Lq ds) (l + ||cr|| L oo (0ii;L9) + ||u|| L oc (0it;i9) ) 



V 



(6.22) 

Plugging (16. 4h into (|6.22j) . we obtain 



^\Ht)\\% 
v 

< 



p 

Then, Young's inequality yields 
fi + A 



t — 

^ (0)Hl 9 + C(p,q,fi)R^J \\a\\ p Lq ds) (1 + \\(r\\ L oo^ t . Lq 



1 FWIILa 



+ gi y (i)jT'lkil£,cfa 



p K 



(6.23) 



for all < t < T. 

Now, we let R be so small that 

C{p,q^)v^R~(l + VR) < l -. 

Due to the fact that ||cr(0)|| W ( K 3) < 5^, we can assume that ||cr(t)||_L<? < \V~R~ in some 
maximal interval (0, i max ) C (0,T). If t max < T, then, ||c(t max )||x9 = \yfR~ and by (|6.23p . 

~VR= \\a(t max )\\ Lq < C(p,q,fi)pR(l + VR) < ^v^R, 
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which is a contradiction. Hence, t max = T and 

\\o-\\lo < ^VE, for all te[0,T]. (6.24) 

Thus, by (16.23p . one obtains (|6.19p . The proof of Lemma 16.31 is complete. □ 
We remark that, from (|6.6p and (|6.19p . one has 

l|V£||£i>(o,T;W(R3)) < C{p,q,fi)R. (6.25) 



7. Global Existence 
In this section, we prove the global existence in Theorem 13.21 Define 

T max := sup jV > : 3 u G W(0,T) with u = Tt(u), such that, 1 1 u| | w(o,T) < R> 
||<S(u) - IHl^ctw 1 .-?) < v 7 ^, lkllz,°°(o,T ; z,9) < V^R, and 

where R was constructed in the previous section. 

If T max = oo, we are done. From now on, we assume that T max < oo. 

7.1. Uniform estimates in time. We now establish some estimates which are uniform 
in time T. First we prove the following energy estimates: 

Lemma 7.1. Under the same assumptions as Theorem \3.1[ we have 

II Vu IIl 2 (o,t ; l 2 (r3)) < CR 2 , (7.1) 

II u IIl°°(0,T;L 2 (]R3)) < CR 2 , (7.2) 

\\ e \\l°°(o,T;L 2 (r 3 )) < CR 2 , (7.3) 

\\P - 1||l°°(o,t ; z, 2 (ir3)) < CR 2 , (7.4) 
where C is a constant independent ofT G (0, r max ). 

Proof. First we recall that 

u G W 1,2 (0, T; L 2 (R 3 )) n L 2 (0, T; W 2 ' 2 (R 3 )) 

and 

p, £ G W 1 ' 2 ^, T; L 2 (R 3 )) n L 2 (0, T; W 1 ' 2 ^)). 

Multiplying equation (|l.lb|) by u, and integrating over M 3 , we obtain, using the conserva- 
tion of mass (jl.lap , 

{\ P H 2 + -^-:(p 1 + l-l)]dx+ I HVu| 2 dx 
dt Jr3 V 2 7 - 1 / Jr3 

pFF T : Vudx. 
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Here, the notation A : B means the dot product between two matrices. Thus, we 

/ (\p\u\ 2 + -^—{p 1 + ~f - l)\ dx + f I p\Vu\ 2 dxds 
Jr3 \2 7 - 1 / Jo Jm* 

= / ( ^Po|u | 2 + —^—r(Po+7- 1)) dx- / / pFF T -.Vudxds. 
From the conservation of mass (jl.laj) , one has 

/ (lp\u\ 2 + —^— -(^-7/9 + 7-1)^) dx+ [ I p\Vu\ 2 dxds 
Jm3\2 7-1 J J J R3 

= [ (^Po|u | 2 + ^-(p2- 7/00 +7-1)) dx- I I pFF T : Vudxds. 
JR3 \2 7-I J Jq J R 3 

On the other hand, due to equations (|l.lc|) and (jl.lap . we have 

= *rr\ F \ 2 + 2 P F '■ ( VuF - u • VF ) 

ot 

= ^|F| 2 + 2/cF : (VuF) - pu • V|F| 2 

= ^|F| 2 + 2/5F : (VuF) + div(/m)|F| 2 - div(pu|F| 2 ) 
= 2pF : (VuF) - div(/9u|F| 2 ). 
Integrating (\7.6h over M 3 , we arrive at 

/ p\F\ 2 dx = pF : (VuF)dx. 

Jm.3 Jrs 



ld_ 

2dt 



Since 



[ [ pF : (VuF)dxds = [ [ pFF T : Vudxds, 
Jo Jm. 3 Jo Jm. 3 

we finally obtain, by summing (|7.5p and (|7.6I) , 

/ f^|u| 2 + ^p|F| 2 + — (p 1 -1P + 1-1)) dx+ [ [ p\Vu\ 2 dxds 
Jus \2 2 7-I J Jo Jrs 

= [ (^o|u | 2 + l-po\Fo\ 2 + —^—(Po-'7Po + J-l))dx. 
Jrs \2 2 7-1 J 

Thanks to Remark 13.41 we have 



p{I + E T ) : Vu = 0. 
Hence, from (jl.lcj) and (jl.lap . we have 

d t (ptrE) = 0. 
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Therefore, from (|7.8p . (|7.9p and the conservation of mass (|l.laj) . we finally arrive at 



/ ( l -p\u\ 2 + \p\E\ 2 + ^-{p< - 1 p + 1 -l)\dx+ I [ p\Vu\ 

JR3 V 2 2 1 - 1 J Jo JR 3 

= I (lpo\^o\ 2 + lpo\Eo\ 2 + ^—(pl-7Po + 7-^)dx<R A . 
JR2 \ 2 2 7 - 1 / 



2 dxds 



(7.10) 



Since p > is a constant and p £ [f , f], then inequalities (|7.ip - (|7.3p follow from (|7.10p . 
and inequality (|7.4p follows from (|7.10p and the following straightforward inequalities: for 
some rj > 0, we have 



x 1 — 1 — 7(x — 1) > 



n\x-l\ 2 , if 7 > 2, 

— 1| 2 , if |x| < 2 and 7 < 2. 

The proof of Lemma 17.11 is complete. □ 

Based on the uniform estimates from Section we have 

Lemma 7.2. Under the same assumptions as Theorem \3.Sk 

ll 5 ( u ) - ^\\l°°(0,T;W 1 ^) < V^R, ||o"||i°°(0,T;i9) < V^R, (7-11) 

/or any T £ [0,T max ]. 

Proof. According to (|6.24p . it is obvious to see that 



max ||erl|L9(i) < VR. 
te[o,T] 



Hence, we are only left to show 



max ||5(u) — 111 (t) < V~R. 

t€[0,T] 

Indeed, for any t £ (0, T), we have, by using (jl.lap and (|6.19p . 
||5(u)(*)-l||2. 



\Po 
\po 



l\\h+J^ s \\S(u)(s)-l\\t q ds 

1112. +7 f (\\s(u)( S )-i\r L ^ 

Q J o 



x 

'■t 



[ \S(u)(s) - l| 9 - 2 (S(u)(s) - l)d s S{u)(s)dx)ds (7 - 12) 

< \\po - n% + - / \\S(u)(s) - lll^ll^^llwds 
1 Jo 

n ( rt (5g-6) P \ E ir / ft \p 

~ q\l " 5(U)(S) ~ 1 " i ' 9 ' 6 dS ) (Jo l|9s5(u)ll "^ S J ' 



where 

(5g-6)(p-l) 

a = h 1. 

3g-6 



(7.13) 
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From (jl.lap and (jl.ldp . we obtain 

II^PlliP(0,T;L9(R 3 )) = l|Vp • u||xp(0,T;L«(K 3 )) 

< 2||(T|| iO o( 0iT;i q( R 3))||u|| L p( 0iT;L oo( K 3)) < C 

On the other hand, from the Gagliardo-Nirenberg inequality, we have 

Hp - i|U, (R3) < c||p - i||i 2(K3) ||v(p - 1)^) < c||p - if L2(R 3) 

with = g|ig. Thus, by Holder's inequality, ((73)1 . and (lfU9D . one has 

" P_1 " L I %^(0 ) T;W(^ - <7 " P_ 1 lli oo (0,T;La(R3))lklfc( ,r;L*(i»s^ ^ 



MIL 



which, together (j7TT2|) and ([7TT3]) . yields 

||«S(u)(i) - 1|U« < Ci2. 

Hence, according to (|6.19p . we obtain, by letting R be sufficiently small, 

max max {\\S(v)(t) - l\\ w i, q(JSL 3 } , ||cr(i)|| L9(B . 3) } < VR. (7.14) 
f(z [o, i \ 

The proof of Lemma 17.21 is complete. □ 
Lemma 7.3. For each 1 < I < 3, ^ satisfies 

d t — + u • V— = - — • VE + V I — — ) E + Vu— + V— (7.15) 



' dxi dxi dxi \ dx\ ) dxi dx 

in the sense of distributions, that is, for all if) £ Cq°(Qt), we have 

dE ^ , , , f T f , ,,dE 
-Otipdxdt + / / div(u^)- 



dxi J J R 3 dxi 

du „ ^ „ / <9u \ ^ „ dE „ du , , , 
- — • V.E + V I - — I E + Vu— + V— ipdxdt, 
o j M 3 v dxi \oxiJ dxi dxi. 



- max ) 



for any T £ (0, T n 

Proof. The proof is a direct application of the regularization. Indeed, one easily obtains, 
using (fL3l) . 

d t (S £ £) + u • V(S £ E) = S e (d t E + u -VE)+u- V(S £ E) - S £ (u ■ VE) 
= S £ (VuE + Vu) + u • V(S £ E) - S £ {u • VE). 
Differentiate (|7.16p with respect to x/, we get 



(7.16) 



ctez J \ dxi 

Notice that 

d / „,^_ N „ , „„A du „ „ / <9u 

<9z; 



(7.17) 



- (u • V(^) - S £ (u • VE)) = ^ • VS £ E -Sj^. VE 
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According to Lemma 6.7 in [26] (cf. Lemma 2.3 in [21]), we know that 



ax; V ox i 



and 



u .V5«,|f)-5«(u.v||)^„, 



in L^O,!"; L q (R 3 )) as e -> 0. Hence, 

c9 



v u-V(S' £J B)-5 E (u-V^) ^0 

dxi 



in L^O, T; L"(M 3 )). Thus, letting e -► in (I7TT71) . we deduce 

d t — + u • V— = -— • VS + V — )E + Vu— + V — , 

OX\ OX\ OX\ \ ox i J OX\ OX\ 

in the sense of weak solutions. The proof of Lemma 17.31 is complete. □ 



Using (|7.15p . formally we have, 



dE 



q-2 



dE , 

dxi 



dxi 

f ( „dE du „ / du\ „ „ dE „du\ 

= / -u-V— - — -VS + V — S + Vu— + V— 

l° o (Qt)II u IIh/ 2 ^II v - e, IIL( 1 ir3) + ll u l 



< CdlVulUoollV^H^^aj + 



dE 



dxi 



q-2 



dE , 

dxi 



< c(l|Vu|| L oc||Vu||^ {R3) + \/^||u|| iy 2, q ||VS||^ ( 1 K3) + 



W 2, q \\X7 EWI^ 



\u\\ W 2,q\\V EW 9 ^ 1 



(7.18) 



We remark that the rigorous argument for the above estimate involves a tedious regular- 
ization procedure as in DiPerna-Lions [9], thus we omit the details and refer the reader to 
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[Sj. Using (jTTLSj) . one obtains 

p 



L'l 



8E(0) 



dxi 
dE{0) 



o ds 



dE 



ds 



L'l 



< 



dxi 
dE{0) p 



+ P - 

Li 1 Jo 



dx. 



+ C\ P - 



L'l 



dE 


p-q 


f 


dxi 


Li 


J» 3 




dE 


p-q 


f 

Jo 


dxi 


Li 



dE 



q-2 



dxi 
Vu 



dE 



+ \/i?||u||l^2, g || V£/||* s /jj,3^ + ||u||v^2,9 IIV-E 1 )!^ 1 



dxi 
ds 



dE 
dxi 



(s) ] dx 



ds 



< 



dE(0) 



dxi 



+ C 



L'l 



Iv^n^tiiVuiiLooiiv^iu, 



+ (1 + Vfl)||u 



\W 2 -i 



Ids 



p-i 
v 



^ +0 (f)U l|V£| H' u 



max IIV-EII + VR+ 1 
te[o,T] 



<5 2p + C ( - ) R p ( max \\VE\\ + Vr+1\. 

te[o,T] 



(7.19) 

Taking the summation over I in (|7.19p and taking the maximum over the time t, one has, 

max \\VE\\ P < 5 2p + C ( - | R p ( max II VE\\ +VR + i) , 
te[0,T] \qj V*e[0,T] / 

and hence, by letting R, 5 be sufficiently small and using (|6.25p . we obtain, 

max \\VE\\ P < 5 2p + CR P < (VR) P . (7.20) 

t€[0,T] 

We are now left to deal with the quantity ||S||^9(r3). To this end, from the Gagliardo- 
Nirenberg inequality, we have 



\E\ 



L'l 



< C\\E\\ e L2 



WE\ 



Li 



with 9 = J^U. Thus, by Holder's inequality, ([73]), and ([OoD 



5q— 6 

l-^ll (5q-6)p 

L zi-s (0,T;L1 



< C\\E\ 



NEW 1 , 



LP(0,T;L<- 



lL°°(0,T;L 2 

Hence, we have the following estimate: 

Lemma 7.4. Under the same assumptions as Theorem \3.Sl it holds 

\\ E \\l^(0,T;Li(R 3 )) < V^) 

for any T E [0,T max ]. 

Proof. By flT3|), QQgp and (f730l) . and letting 

(5g — 6)(p — 1 



(7.21) 



(7.22) 



O: 



3q-6 



+ 1, 
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one obtains, 
\\E(t)\\h 

'* d 





+ f 
Jo 


= \\E(0)\\% 


a 
+ - 

q 


= \\E(0)\\h 


a 
+ - 

q 


< 


a 
+ - 

q 


< \\E(0)\\% 


a 
+ - 

q 



ds 



\E(s)\\t q ds 



J (\\E{s)\\ a ~ q j \E(s)\ q - 2 E{s) VuE + Vu-u-VE dx\ ds 



2||Vu||loo \\E\\l9 + ||Vui|£9 



p-1 

t (5qr-6)p \ p 

\\E(s)\\ L p- 6 dt) ||u|| LP(0iT;H/2 „ 



ds 

(7.23) 



x 2 sup \\E(t)\\ Lqm + 1 
V te(o,T max ) 



<\\e(o)\\% + [y) r [J l|s(s)ll ^~ 6 dt ) 



< ||£(0)||g, +CR\\E\\%\_ 6)p 

L a?-B (0,T;L9 

Then, according to (|7.21[) . one has, for all i G [0, T max ], 

w < ^ + CR a < yfit, (7.24) 

if R is sufficiently small. Thus, (|7.22p follows from (|7.24p . The proof of Lemma 17.41 is 
complete. □ 

Lemma 17.41 together with (|7.20p and Lemma 17.21 gives 

max max{||cS(u) - l|| w i,,(t), ||o-|k,(*), llTfu)!!™-!,, (t)\ <CR<\^R. (7.25) 

te[o,T] L ' 

Similarly, we can obtain 

max max (||<S(u) — l|| w i,2(i), (t), \\T(u)\\ wl , 2 (t)} <CR<VR. (7.26) 

te[o,T] 

7.2. Refined estimates on Vp and V-E. In order to prove Theorem 13.21 we need some 
refined estimates on || V/0[|£2( O) t;L9(K 3 )) an d \\^E\\ L 2^^ T . Lq ^ 3 )y 

Lemma 7.5. 

l|Vp|U 2(0 ,T;L 9 (M3)) < R 2 , (7.27) 

for any T £ (0,T max ). 

Proof. Taking the divergence in (ll.lbj) . and using divu = 0, one obtains 

AP(p) = div(div(p££ T )) + div(div(p£)) - div(pu • Vu) - div((p - l)d t u). (7.28) 
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Since, div(/?(I + E) T ) = 0, we get 



div(div(p£)) = — — {p Ei3 ) = { pEi3 ) = -Ap 



d d ^ d d 
dxi dxj 3 dxj dxi 
in the sense of distributions. Hence, (|7.28[) becomes 

AP(p) + Ap = divdiv( P EE T ) - div(pu • Vu) - div((p - l)d t u). (7.29) 
Hence, one obtains, using L q theory of elliptic equations and Taylor's formula, 

II v pIIl 2 (o,t ; l<?(r3)) < c(llP u Vu|| L 2( 0iT;L9(i >3)) + ||div(p£ , £ ,T )|| L 2( 0)T . L3 ( K 3)) 

+ || (p - l)<9tU|| i 2( 0iT;L g( R 3)) 

< c(||/ 5 IIl° o (Qt)II^ U IIl 2 (0,T;L°°(M 3 ))II u IIl°°(0,T;L9( 

+ I|Vp||l2( jyl9(r3))||£|||oo ( q t) 

+ \\p\\l^(Qt)\\^ E \\l 2 (0,T;Li(R3))\\E\\ L oo^q t ^ 
+ \\P— 1 IIl°°(Qt)H < ^ u IIl 2 (0,T;L<J 

< C(l|p||L°°(Q T )ll u llL 2 (0,T;W 2 >9(R3))||u|| L oo( 0iT;L q( R 3)) 

+ H^llL 2 (0,T;i9(R 3 ))ll^'lll, oo (0,T;W 1 .«(R 3 )) 

+ IIpIU°°(Q t )||V£?||£2(o ) T;L9(R3))||^|U'»(0,T;W 1 .«(I11 3 ))) + ^ 

< CR 2 . 

The proof of Lemma 17.51 is complete. □ 

Ir order to refine [|V-E[|£2(o,t;.m(K 3 ))j we need the following estimate: 
Lemma 7.6. 

3 

\\dtu\\ L 2 {0)T . Lq{R3)) < Cm, (7.30) 

for any T G (0,T max ). 

Proof. We first notice that, by the Gagliardo-Nirenberg inequality, for q £ (3,6], 

II u IIl 2 (0,T;L9(R3)) < C|l Vu llL 2 (O,T;L 2 (R3))ll u lll2( OiT;i 2( R 3)) < CR l+9 , (7.31) 

with 6 = G (3, !]• Next, we multiply (jl.lbp by d 4 u and integrate over R 3 x (0,i) to 



deduce 

r-t 



[ [ p\d t u\ 2 dxds + n I \Vu{t)\ 2 dx 

JO JR3 JR3 

= /W \Vuo\ 2 dx — / / pu • Vu • dtudxds — / / VPdtudxds 

JR3 Jo VR3 JO JR 3 

+ / / div(pEE T )d t udxds + / / div(pE)d t udxds 
Jo Jr 3 Jo Jr 3 

f 4 
:=H / |Vu | 2 dx + VJi, 



(7.32) 
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with the following estimates on 7j (i = 1...4): recalling Qt = R 3 x (0, T), 

\h\ < ll^p5iU|| L 2(Q T )||^/p|| i oc(Q T )||u|| Lo o( 0iT;i 2( R 3))||Vu|| i 2( 0iT;i o 

^ l ( I P\d t u\ 2 dxds + CR 4 ; 
8 Jo Jr 3 

i r T r 

p\d t u\ 2 dxds + C\\Vp\\l 2{QT) 



o Jr 3 



\h\ < C\\V p\\ L 2 {QT) \\^d t u\\ L 2 (QT) < 

^ I I I p\d t u\ 2 dxds + CR 4 ; 
8 Jo Jr 3 

\h\ < C\\Vp\\ L 2 (QT) \\E\\ 2 Lx{QT) \\^d t u\\ L 2 {QT) + C||E|| L ^ ( Q T) ||V J B|| i 2 ( Q T) ||^9iu|| L 2 ( Q T) 



< - / / p\d t u\ 2 dxds + CR A ; 
8 Jo Jr 3 



\h\ < 



JR 3 



d t {pE)Vu.dxds 



+ 



+ 



p(T)E(T)Vu(T)dx 



p E Vu dx 

< (IIpIIl-(q t )II^^||l2 (Qt) + \\E\\ LX{QT) \\d t p\\ L 2 {QT) )\\Vu\\ L 2 {QT) 

+ cr 3 + (\\v P (T)\\ L 2 m \\E(T)\\ Loam + iiv^^ii^^n^r)!!^^)^^)!!^^) 

< CR 3 , 

where, for the estimate 1^, we used equations (|l.la|) . (|l.lcj) . Lemma [7,11 and estimate 
(|7.31|) . Thus, from (|7.32p . one obtains 

\\d t u\\ L 2 (QT) <CRl, (7.33) 

and 

3 

\\^\\L^(0,T;L^(R^ < CR72. (7.34) 

Now, we differentiate (jl.lbj) with respect to t, multiply the resulting equation by dtu, and 
integrate it over M 3 to obtain 

--j- / p\dtu\ 2 dx + p f \Vdtu\ 2 dx 

2 dt J R 3 J R 3 

= — j dtp\dtu\ 2 dx — j dtpu ■ Vu • dtudx 
2 Jr 3 Jk 3 

— / pd t u ■ Vu • d t udx — j pu- Vdtu ■ dtudx — / VdfPdtudx (7 
Jr 3 Jr 3 Jr 3 {(.60) 

- [ d t (pEE T )V d t udx - [ d t (pE)Vd t udx 
Jr 3 Jr 3 

where using (|7.34|) . we can control Jj (j = 1...7) as follows: 



Vpu\dtu\ 2 dx 



< ||9tu||| 6 ||Vp|| L 3||u|| L 3 < <7i? 2 ||Vd t u||| 2 ; 



and 
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\H < ||$u|| iB ||Vp|| L 3||u||2 8 ||Vu|Ue < #||V3 t u|| L2 ||Vu|| 2 2 ||Au|| L2 

< R 2 \\Vd t u\\ 2 L2 +i? 3 ||Au||2 2 ; 

|J 3 | < ||p||l-||Vu|| l 3 ||^u||| 6 < ||p|U»||u||| 2 ||Vu||| 2 ||Va t u||| 2 < CRl\\Vd t u\\l 2 ; 
\Ja\ < ||p||L^||u|| L 3||V^u|| i2 ||^u|| L 6 < CR\\Vd t u\\ 2 L2 ; 
\Js\ < C\\d t p\\ L 2\\Vd t u\\ L 2 < C\\Vpu\\ L2 \\Vd t u\\ L 2 

< C\\Vp\\ L s\\n\\ L e\\Vd t u\\ L 2 < C||Vp|| L 3||Vu|U 2 ||va t u|U 2 

<|||V^u|l! 2 + Ci? 2 ||Vu||| 2 ; 

|j 6 | < WdtphzWEf^wvdtuWv + \\p\\ L ~\\E\\ L oo\\d t E\\ L 4vd t u\\ L2 

< R 2 \\Vp\\ L 3 ||u|| L6 ||Va t u|| L2 + R\\ VuE -U-X7E + Vu|| L 2 \\Vd t u\\ L 2 

<E 2 ||Vp|| L3 ||Vu|| i2 ||VeHlz3 
+ R(\\Vu\\ L3 \\E\\ L6 + ||V^|Us||u|| z8 + ||Vu|| L2 ) ||Vd t u|| L2 

< fl 6 ||Vu|| 2 2 + ^\\Vd t u\\l 2 + i? 2 (i* 2 ||Vu||f 3 + ||Vu|| 2 2 ); 

\Jr\ < ||p||l-||^|| l2 ||V^u|| l2 + \\E\\ L ~\\d t p\\ L2 \\Vd t u\\ L2 

< ^\\Vd t u\\l 2 + i? 2 ||Vu|| 2 3 + ||Vu|| 2 2 +R 2 \\Vpu\\l 2 

< ^\\Vd t u\\ 2 L2 + i? 2 ||Vu|| 2 i3 + ||Vu|| 2 2 + fl 4 ||Vu|| 2 2 . 

We remark that in the above estimates, we used several times the interpolation inequality: 

II/I|W 2 .3(R3) < ||/||VF 2 , 2 (R3)||/||^ 2 , 9 ( R 3) 

for some 9 £ (0, 1). These estimates and (|7.35p imply that, for R sufficiently small, 

— / p\d t u\ 2 dx + ^ [ |V5 4 u| 2 ^<i? 3 ||Au|| 2 2 +C||Vu|| 2 2 +i? 2 ||Vu|| 2 3 . (7.36) 
^ at Jr3 8 J R 3 

Integrating (|7.36p over (0, T), we obtain that, using (|7.1|) . 

\\Vd t u\\ LHQT) < CRl. (7.37) 

Here we used the estimate 

\\p d t u(0)\\ L 2 < C\\u • Vu || L2 + ||Au || L2 + \\Vp \\ L 2 + \\VE \\ L 2 < 5 4 

by letting t = in (|l.lbj) . Thus, by (|7.33[) . (|7.37p and the Gagliardo-Nirenberg inequality, 
we obtain 

II^u||l2(o,T;M(R3)) < II^ u IIl2(q t )II V ^ u II12 ( q t) < crI, 
for some 6 G (5, 1]. The proof of Lemma 17.61 is complete. □ 

With (|7.30p in hand, we can now get the estimate for \\VE\\ L 2^ 0>T . Lq ^3^. 
Lemma 7.7. 

I|V£|| L 2(o,t ; l*(k 3 )) < CRi , (7.38) 

for any T E (0,T max ). 
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Proof. Substituting the following two facts 

d t divE = div(-u • VE + VuE) + Au, 

and 

dW(pE) = div((p - l)E) + div£, 

into (|l,lb|) . multiplying the resulting equation by |divS| 9_2 div£' and integrating it over 
R 3 , we can obtain 

p d 



q dt 

< 

+ 
+ 
+ 



div^li, + lldiv^Hi, 
pd t u\drvE\ q ~ 2 divEdx 
VP\divE\ q - 2 divEdx 



+ 



+ 



puVu\dwE\ q ~ 2 dwEdx 
I div( P EE T )\divE\ q - 2 divEdx 



div((p- l)E)\dWE\ q - 2 divEdx 



/ div(VuE-u- VE)\divE\ q ~ 2 dwEdx 
Jr 3 



(7.39) 



M ™, 



m=l 



where 



Mi < ||p|M|3tu|| L9 ||div£|||/; 



M2 < ||p||l°° || u I|l9 II Vu||ioo HdivE 1 !!^ 1 < i?||u||^2, 

Ma^CHVpll^lldiv^ll^ 1 ; 
M 4 < ||Vp|| L9 ||£||! 00 ||div£||!~ 1 + || / o|| L oo||E|| L oo||V J B|| L9 ||div^||^ 1 
< (R 2 \\Vp\\ Lq + R\\VE\\ Lq )\\divE\\ q - 1 ; 
M 5 < ||p- l|| L oo||V£|| L g||div£||'^ 1 + ||V / o|| L g||E|| L oc||div£'||^ 1 



Idiv^Hl; 1 ; 



<i2(||V£|| L9 + ||Vp|| L 0l|div£|| 9_1 - 



\Li ' 



M 6 < (||Vu|| LO o||V^||l9 + [|Au[|r«||^|U'»)||divS||2 i , 1 < i?||u 



W 2 'i 



\d\vE 



19-1 



Here in Mq, we used divu = 0. With those estimates in hand, we multiply (|7.39p by 
|div.E| L ~ 9 to deduce that, using Young's inequality, 

^j t \\divE\\ 2 Lq + \\divE\\ 2 Lq < C\\d t u\\ 2 Lq + R 2 \\u\\ 2 w2 , q + || Vp||| 9 + R 2 \\VE\\ 2 Lq . (7.40) 

On the other hand, we still have 

||curLE|||, < \\E\\ 2 L oo\\VE\\ 2 Lq < CR 2 \\VE\\ 2 Lq . 
Hence, substituting this into (|7.40p . we get 

l^lldiv^Hl, + HV^Hi, < C\\d t u\\ 2 Lq + R 2 \\u\\ 2 w2 , q + ||Vp||| 9 + CR 2 \\VE\\ 2 Lq . (7.41) 
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Integrating (|7.4ip over (0, T) and using the estimates (|7.27p . (|7.30p . we obtain 

3 



\L 2 (0,T;L1( 

The proof of Lemma 17.71 is complete. □ 
7.3. Proof of Theorem I3.2L First we have the following estimate on ||u|| w(o,T) : 
Lemma 7.8. 

l|u||w(o,r) < R, (7.42) 

for all T G [0 

i?max]) if R is small enough. 
Proof. First from Theorem 14.11 we have 

ll u llw(0,T) = ll^( u )llw(0,T) 

< C(q)(\\u \\ v 2, q + ||(1 - <5(u))a t u|| L 2( 0iT . ig(K 3)) 

+ ll«5(u)(u • V)u|| L 2 (0jT;iq(R 3 )) + ||VP(5(u))|| i2(0iT;i9 
+ ||div(«S(u)(7 + T(u))(7 + T(u)) T )\\ L 2^ T . LHm) ] 

From the previous computations in Lemma 17.21 we have 

||(1 - S(u))dtu|| L 2( 0)T;L ,( K 3)) < ||«5(u) - 1||l°°(q t )II u IIw(o,t) 

< C||<S(ll) - 1||loo(0 ) T;W1.9(R3))||u||w(0,T) 

< C(q)R?. 

Similarly, one has 

||«S(u)(u- V)u|| L 2 (0iT;L9(R 3 )) < ^(u)!!^^)!^!!^^ < C(q)R 2 ; 

\\VP(S(u))\\ L 2 (mm3)) = \\P'(S(u))Vp\\ L 2 ( , T . Lm3)) < CR 2 ; 



and 

||div(«S(u)(7 + T(u))(7 + T(u)) T )|| L2(0iT;L9 

< \W P\\l 2 {0,T;L1(M?)) (l + ll-^ll £°°(Q T )) 

+ W E \\l 2 (0,T;L1(TS. 3 ))\\p\\l x (Qt) (l + \\E\\ L oo(Q T ^) 

< C(q)R*. 

Thus, we obtain, 

||u|| w( o,t) < C(q)(R 2 +Rl)<R, for all T G [0, T max ], 

if R is small enough. The proof of Lemma 17.81 is complete. □ 

Finally we are in the position to give the proof of Theorem 13.21 as follows. 
Let T n / 

Tmax be an increasing sequence with limit equal to T" max . Then, it is necessary 
to have, according to Theorem 14.21 

sup ||u(T n )|| x <£(dR, 
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which, combining with the definition of T max , implies 

\\(u(T n ),S(u)(T n ) - l,T(u)(T n ))\\ vr < C(q,c )R. 

Therefore, if we take each triple (u(T n ), <S(u)(T n ) — 1, T(u)(T n )) as a new initial condition, 
then according to Theorem 13, 1\ we know the solution can be extended to the interval 
(0, T„+To), for some To > 0. Hence, the solution can be extended to the interval (0, T max + 
T ), which, together with (TT25]) and (f7726l) . is a contradiction with the maximality of T max . 
So it cannot be true that T max < oo, but T max = oo, i.e. u, 5(u), and T(u) are well-defined 
on M 3 x (0, oo). The proof of Theorem 13.21 is complete. 

Acknowledgments 

Xianpeng Hu's research was supported in part by the National Science Foundation grant 
DMS-0604362 and by the Mellon Predoctoral Fellowship of the University of Pittsburgh. 
Dehua Wang's research was supported in part by the National Science Foundation under 
Grant DMS-0604362, and by the Office of Naval Research under Grant N00014-07-1-0668. 

References 

[1] Admas, R. A., Sobolev spaces. Pure and Applied Mathematics, Vol. 65. Academic Press, New York- 
London, 1975. 

[2] Amann, H., Linear and quasilinear parabolic problems. Vol. I. Abstract linear theory. Birkhuser 

Boston, Inc., Boston, 1995. 
[3] Bergh, J.; Lofstrom, J., Interpolation spaces. An introduction. Grundlehren der Mathematischen 

Wissenschaften, Springer- Verlag, Berlin-New York, 1976. 
[4] Bird, R. B.; Armstrong, R. O; Hassager, O., Dynamics of polymeric liquids. Vol. 1. Fluid mechanics. 

Wiley Interscience, New York, 1987. 
[5] Chemin, J.; Masmoudi, N., About lifespan of regular solutions of equations related to mscoelastic 

fluids. SIAM J. Math. Anal. 33 (2001), 84-112. 
[6] Chen, Y.; Zhang, P., The global existence of small solutions to the incompressible viscoelastic fluid 

system in 2 and 3 space dimensions. Comm. Partial Differential Equations 31 (2006), 1793-1810. 
[7] Dafermos, C. M., Hyperbolic conservation laws in continuum physics. Second edition. Grundlehren 

der Mathematischen Wissenschaften, 325. Springer- Verlag, Berlin, 2005. 
[8] Danchin, R., Density-dependent incompressible fluids in bounded domains. J. Math. Fluid Mech. 8 

(2006), 333-381. 

[9] DiPerna, R. J.; Lions, P.-L., Ordinary differential equations, transport theory and Sobolev spaces. 
Invent. Math. 98 (1989), 511-547. 

[10] Goktekin, G. B.; Bargteil, A. W.; O'Brien, J. F., A method for animating viscoelastic fluids, ACM 
Transactions on Graphics 23 (2004), 463-468. 

[11] Gurtin, M. E., An introduction to Continuum Mechanics. Mathematics in Science and Engineer- 
ing, 158. Academic Press, New YorkLondon, 1981. 

[12] Joseph, D., Fluid dynamics of viscoelastic liquids. Applied Mathematical Sciences, 84. Springer- 
Verlag, New York, 1990. 

[13] Kunisch, K.; Marduel, M., Optimal control of non-isothermal viscoelastic fluid flow, J. Non- 
Newtonian Fluid Mechanics 88 (2000), 261-301. 

[14] Larson, R. G., The Structure and Rheology of Complex Fluids. Oxford University Press, New York, 
1995. 

[15] Lei, Z.; Liu, C; Zhou, Y., Global existence for a 2D incompressible viscoelastic model with small 

strain. Commun. Math. Sci. 5 (2007), 595-616. 
[16] Lei, Z.; Liu, C; Zhou, Y., Global solutions for incompressible viscoelastic fluids. Arch. Ration. Mech. 

Anal. 188 (2008), 371-398. 



STRONG SOLUTIONS TO INCOMPRESSIBLE VISCOELASTIC FLUIDS 



41 



[17] Lei, Z.; Liu, C; Zhou, Y., Global solutions for compressible viscoelastic fluids with small initial data. 
Preprint. 

[18] Lei, Z.; Zhou, Y., Global existence of classical solutions for the two-dimensional Oldroyd model via 

the incompressible limit. SIAM J. Math. Anal. 37 (2005), 797-814. 
[19] Lin, F; Liu, C; Zhang, P., On hydrodynamics of viscoelastic fluids. Comm. Pure Appl. Math. 58 

(2005), 1437-1471. 

[20] Lin, F; Zhang, P., On the initial-boundary value problem of the incompressible viscoelastic fluid 

system. Comm. Pure Appl. Math. 61 (2008), 539-558. 
[21] Lions, P. L., Mathematical Topics in Fluid Mechanics. Vol. 1. Incompressible Models. Oxford Lecture 

Series in Mathematics and its Applications, 3. Oxford Science Publications. The Clarendon Press, 

Oxford University Press, New York, 1996. 
[22] Lions, P. L.; Masmoudi, N., Global solutions for some Oldroyd models of non-Newtonian flows. 

Chinese Ann. Math. Ser. B 21 (2000), 131-146. 
[23] Liu, C; Walkington, N. J., An Eulerian description of fluids containing visco-elastic particles. Arch. 

Ration. Mech. Anal. 159 (2001), 229-252. 
[24] Matsumura, A.; Nishida, T., The initial-value problem for the equations of motion of viscous and 

heat-conductive gases. J. Math. Kyoto Univ. 20 (1980), 67-104. 
[25] Matsumura, A.; Nishida, T., Initial-boundary value problems for the equations of motion of com- 
pressible viscous and heat- conductive fluids. Comm. Math. Phys. 89 (1983), 445-464. 
[26] Novotny, A.; Straskraba, I., Introduction to the mathematical theory of compressible flow. Oxford 

Lecture Series in Mathematics and its Applications, 27. Oxford University Press, Oxford, 2004. 
[27] Oldroyd, J. C, On the formation of rheological equations of state, Proc. Roy. Soc. London, Series A 

200 (1950), 523-541. 

[28] Oldroyd, J. G., Non-Newtonian effects in steady motion of some idealized elastico-viscous liquids, 

Proc. Roy. Soc. London, Series A 245 (1958), 278-297. 
[29] Renardy, M.; Hrusa, W. J.; Nohel, J. A., Mathematical Problems in Viscoelasticity. Longman Scientic 

and Technical; copublished in the US with John Wiley, New York, 1987. 
[30] Salvi, R.; Straskraba, I., Global existence for viscous compressible fluids and their behavior as t — > oo. 

J. Fac. Sci. Univ. Tokyo Sect. IA Math. 40 (1993), 17-51. 
[31] Sideris, T. O, Nonlinear hyperbolic systems and elastodynamics. Phase space analysis of partial 

differential equations. Vol. II, 451-485, Pubbl. Cent. Ric. Mat. Ennio Giorgi, Scuola Norm. Sup., 

Pisa, 2004. 

[32] Sideris, T. O, Personal notes on viscoelasticity, 2007. 

[33] Sideris, T. O; Thomases, B., Global existence for three-dimensional incompressible isotropic elasto- 
dynamics via the incompressible limit. Comm. Pure Appl. Math. 58 (2005), 750-788. 

[34] Slemrod, M., Constitutive relations for Rivlin-Ericksen uids based on generalized rational ap- proxi- 
mation. Arch. Ration. Mech. Anal. 146 (1999), 73-93. 

[35] Yu, J. D.; Sakai, S.; Sethian, J. A., Two-phase viscoelastic jetting , J. Computational Physics 220 
(2007), 568-585. 

Department of Mathematics, University of Pittsburgh, Pittsburgh, PA 15260, USA. 
E-mail address: xihl50pitt.edu 

Department of Mathematics, University of Pittsburgh, Pittsburgh, PA 15260, USA. 
E-mail address: dwang@math.pitt.edu 



